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LIST  OF  ABBREVIATIONS 


Av  “  surface  area  of  water  droplets  per  voliune  of  cooling  tower 

Cp,  -  constant  pressure  specific  heat  of  moist  air 

Cpw  -  constant  pressure  specific  heat  of  liquid  water 

C,  -  derivative  of  saturation  air  enthalpy  with  respect  to 

temperature 

h,  -  enthalpy  of  moist  air  per  mass  of  dry  air 

he  -  convection  heat  transfer  coefficient 

hp  “  diffusion  mass  transfer  coefficient 

hg  -  enthalpy  of  water  above  reference  state  for  liquid  water 

at 

h,  -  enthalpy  of  saturated  air 

m*  -  mass  flow  rate  of  water 

m,  -  mass  flow  rate  of  air 

m*  -  ratio  of  air  flow  rate  to  water  effective  capacitance  rat3 
for  effectiveness  analysis 

Le  -  Lewis  number 

Ntu  -  overall  number  of  transfer  units 

Q  -  actual  heat  transfer  rate 

T„f  -  reference  temperature  for  zero  enthalpy  of  liquid  water 
(32  °F) 

T„  -  water  temperature 

V  -  volume 

e,  -  air-side  heat  transfer  effectiveness 

w,  -  air  humidity  ratio 

w,  -  humidity  ratio  of  saturated  air 
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a  -  air  stream  conditions 

i  -  inlet  state  conditions 

o  -  outlet  state  conditions 

T  -  total 

w  -  water  stream  conditions 
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Optimization  of  the  design  of  chiller/cooling  tower  systems 
for  large  office  buildings  require  an  accurate  detailed  model  for 
the  cooling  tower.  This  model  must  be  precise  and  allow  for 
varying  air  and  water  flow  rates,  ambient  air  conditions  and 


entering  water  temperatures.  The  optimization  computational 
requirement  is  very  time  consuming,  however,  due  to  the  tower  model 


requiring  solution  of  three  simultaneous  differential  equations. 

The  purpose  of  this  project  was-t^  develop^  a  computer  program 
utilizing  an  effectiveness  model  to  reduce  computational  time  while 
maintaining  good  correlation  to  existing  programs.  The  methodology 
involved  the  use  of  a  tower  model  program _d^eloped  by  Eric  Weber 
[1988]  and  the  modifications  made  unto  by  Thane  Joyce  [1990]  that 
integrated  a  centrifugal  chiller  model,  cooling  tower  model,  and 


system  optimization  strategy.  The  effectiveness  model  generated  by 
James  Braun  [1988]'was  utilized  in  place  of -Weber^^ tower  model. 

The  results  are  very  encouraging.  The  computational  time 
required  to  optimize  the  design  was  reduce  approximately  75%  while 
maintaining  an  accuracy  of  tower  water  temperature  range  of 
approximately  15%.  This  reduction  in  computer  time  will  allow  an 
engineer  the  flexibility  to  examine  many  systems  and  aetermine  the 
optimum  design  in  a  shorter  time  frame. 
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INTRODUCTION 


The  design  of  heating,  ventilation  and  air  conditioning  (HVAC) 
for  large  buildings  tasks  the  engineer  with  manipulating  a  plethora 
of  design  parameters.  With  limitations  on  funding  and  typically  a 
short  fused  deadline,  engineers  rely  on  rule-of-thumb  or  past 
experience  when  designing  these  systems.  No  attempt  is  made  to 
optimize  the  system  as  a  unit  bringing  together  all  component 
parts:  pumps,  chiller,  cooling  tower,  pipe  size,  etc. 

An  optimization  methodology  for  the  chiller/cooling  tower 
system  was  developed  by  Thane  Joyce.  The  results  of  his  work  show 
that  an  optimized  system  can  save  both  construction  and  operating 
costs.  The  computational  requirements  to  perform  the  optimization, 
however,  is  very  time  consuming  and  therefore  the  methodology  was 
not  applicable  as  a  design  tool.  In  order  to  make  the  program  more 
useful,  the  computational  time  needed  to  be  reduced.  This  is 
accomplished  by  employing  Braun's  effectiveness  model  in  lieu  of 
the  Weber  cooling  model  in  Joyce's  program.  The  resulting  program 
was  then  compared  to  the  original  Joyce  program  for  accuracy  and 
reduction  in  computational  time. 
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DEVELOPMENT  OF  EFFECTIVENESS  COMPUTER  MODEL 


The  effectiveness  model  used  to  replace  the  tower  model 
developed  by  Weber  is  based  on  the  work  by  Jeunes  Braun.  A 
schematic  of  a  counterflow  cooling  tower  is  provided  in  figure  1. 


Figure  1,  Schematic  of  a  Counterflow  Cooling  Tower  [Braun, 
1988] 
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The  following  assumptions  were  made  by  Braun  to  derive  the  basic 
model  equations: 

1.  Heat  and  mass  transfer  in  the  direction  normal  to  flows  only. 

2.  Negligible  heat  and  mass  transfer  through  tower  walls  to 
environment. 

3.  Negligible  heat  transfer  from  the  tower  fans  to  the  air  or 
water  streams. 

4.  Constant  water  and  dry  air  specific  heats. 

5.  The  mass  fraction  of  water  vapor  in  the  mixture  of  air  and 
vapor  is  approximately  equal  to  the  humidity  ratio. 

6.  Uniform  temperature  throughout  the  water  stream  at  any  cross 
section. 

7.  Uniform  cross-sectional  area  of  the  tower. 

From  steady-state  energy  and  mass  balances  on  an  incremental 
volume,  refer  to  figure  1,  the  following  differential  equations  are 
developed  by  Braun  and  Weber: 
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Knowing  the  inlet  conditions,  the  number  of  transfer  units,  and  the 
Lewis  number,  the  exit  state  of  both  the  air  and  water  streams  can 
be  determined.  Using  these  same  equations,  Weber  took  an  iterative 
approach  with  respect  to  the  outlet  humidity  ratio  and  the  exit 
water  temperature,  numerically  integrating  the  equations  over  the 
entire  volume  from  air  inlet  to  outlet. 

The  effectiveness  model  simplifies  this  analysis  considerably. 
By  employing  Merkel's  assumptions  of  neglecting  water  loss  and 
setting  the  Lewis  number  equal  to  unity,  the  three  differential 
equations  reduce  to; 


dV 


(4) 


^  m„(dhJdV) 


Equation  (5)  can  be  written  in  terms  of  air  enthalpies  only  by 

introducing  a  new  variable,  C,,  termed  the  saturation  specific  heat: 

dV  KC,{dhJdV) 

~dr  ■  — wc - 
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where 


dh 

“•*  r-r. 


C,  is  equivalent  to  the  derivative  of  the  saturated  air  enthalpy 
with  respect  to  temperature  evaluated  at  the  water  temperature  and 
has  the  units  of  specific  heat.  Braun  indicates  that  if  the 
saturation  air  enthalpy  were  linear  with  respect  to  temperature 
then  the  exit  conditions  could  be  solved  for  analytically.  His 
figure  2,  shown  below,  depicts  the  variation  of  the  saturation 
enthalpy  with  temperature  for  typical  water  inlet  and  outlet  points 
along  with  a  straight  line  connecting  the  points. 


Figure  2.  Saturation  Air  Enthalpy  versus  Temperature  [Braun, 
1988] 
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It  is  clear  form  the  figure  that  the  saturation  air  enthalpy  does 
not  vary  linearly  with  temperature,  however,  by  choosing  an 
appropriate  average  slope  between  the  inlet  and  outlet  states,  an 
air-side  effectiveness  relationship  in  terms  of  C,  can  be  developed. 
This  air-side  effectiveness  is  defined  as  the  ratio  of  the  actual 
heat  transfer  to  the  maximum  possible  heat  transfer  if  the  exiting 
air  were  saturated  at  the  temperature  of  the  incoming  water  (i.e., 
^•,0  ~  • 

Writing  the  heat  transfer  in  terms  of  this  effectiveness 
yields; 

6  =  <»> 


where 


€g  -  air-side  heat  transfer  effectiveness 
1  -  exp  (  -  Ntu  (1  -  m  * ) ) 

1  -  m  *  exp  (  -  Ntu  (1  -  m*)) 


(9) 


and 


m  = 


m. 


(10) 


The  average  value  for  the  saturation  specific  heat  is  estimated  as 
the  average  slope  between  the  inlet  and  outlet  state  points: 


~T 
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(11) 
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The  outlet  air  enthalpy  and  water  temperature  can  then  be 
calculated  using  the  following  equations: 


^a,o  ^aj  ~ 


(12) 


’  w,o 


=  <<  -  K  iKo  -  Ki)  +  r 


m  C 
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(13) 


In  the  computer  model  employing  the  air-side  effectiveness,  the 
water  loss  is  neglected,  so  m^j  is  set  equal  to  m^^.  Since  C, 
depends  upon  the  solution  for  the  exit  conditions  is  iterative 
r-ing  equations  9  through  13.  The  saturation  specific  heat,  C,,  is 
only  weakly  dependent  upon  T^^,  however,  so  any  reasonable  initial 
guess  for  T^,,  (such  as  the  inlet  wet  bulb  temperature  of  the  air) 
typically  results  in  convergence  in  only  two  iterations. 
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CONCLUSIONS 


The  results  of  this  exercise  are  very  promising.  The 
computational  time  savings  for  the  effectiveness  method  was  as  much 
as  83%  over  the  Joyce  method,  while  the  percent  error  in 
calculating  the  energy  consumed  was  less  than  14%  for  all  cases. 

The  first  portion  of  the  project  is  dedicated  to  evaluating 
the  affect  of  employing  the  effectiveness  model  on  the  delta 
temperature  across  a  counterflow  cooling  tower.  Weber's  computer 
program  is  used  to  generate  the  control  data  set.  The  revised 
program  generated  the  test  data  set.  Figure  3  displays  the  results 
of  this  evaluation.  Numerical  data  is  found  in  Appendix  B. _ 

ACCURACY  OF  WEBER  vs  EFFECTIVENESS 

INCLUDING  +  &  -  15%  CONFIDENCE  BANDS 
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Figure  3.  Accuracy  of  Weber  versus  Effectiveness  Model 
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The  effectiveness  model  computer  program  is  provided  as  Appendix  A. 
This  model  gives  a  larger  delta  temperature  across  the  cooling 
tower  in  all  cases  compared  with  the  exact  data  set.  One  reason 
why  the  effectiveness  model  overpredicts  the  tower  delta 
temperature  is  that  Weber's  model  is  developed  :for  a  crossflow 
cooling  tower  while  Braun's  analysis  is  for  a  counterflow  cooling 
tower.  If  improved  accuracy  in  calculating  the  delta  temperature 
across  the  tower  is  desired,  a  crossflow  correction  factor  can  be 
applied  to  the  effectiveness  model.  Linearization  of  saturation 
enthalpy  between  the  inlet  and  outlet  points  tends  to  overpredict 
the  tower  heat  transfer  and  is  another  factor  contributing  to  this 
error . 

To  test  and  compare  the  effectiveness  model  to  the  Weber 
cooling  tower  model  in  Joyce's  program,  several  variables  are 
manipulated  while  holding  others  constant.  The  variables 
investigated  are:  1)  air  flow  rate  through  tower  (cfm) ,  2)  tower 
size,  and  3)  pump  size  (gpm) .  These  variables  are  considered  the 
primary  parameters  for  designing  cooling  towers.  Figures  4  through 
9  show  the  results  of  the  comparisons  for  these  various  parameters. 
The  results  are  encouraging  and  in  each  case  the  trends  for  both 
models  are  similar.  Appendix  C  contains  the  numerical  data  used  in 
generating  these  figures. 

Figure  4  depicts  the  plot  of  effectiveness  data  versus  Joyce 
data  for  various  air  flow  rates  through  the  cooling  tower  (CFM) . 
The  trends  in  both  computational  methods  for  this  design  parameter 
are  very  similar. 
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ACCURACY  OF  JOYCE  VS  EFFECTIVENESS 

FOR  VARIOUS  CFWS 
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The  time  savings  for  this  variable  are  substantial,  as  shown  in 
figure  5. 

Six  different  tower  air  flow  rates  were  used  while  holding  the 
other  variables  constant.  The  maximum  error,  of  12.8%,  occurs  at 
a  flow  rate  of  160000  CFM.  The  time  savings  are  between  33  and  37 
minutes,  or  roughly  81%  improvement  over  the  more  exact  Joyce 
model. 

The  next  variable  manipulated  is  tower  size.  Actual 
manufacture's  performance  data  for  Baltimore  Aircoil  series  3000  is 
used  as  input  parameters.  The  results  are  very  favorable  giving 
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Figure  5.  Time  Savings  at  Various  CFMs. 


similar  trends  for  both  models,  figure  6.  The  time  savings  using 
the  effectiveness  model  over  the  Joyce  method  for  this  parameter 
are  also  quite  impressive,  refer  to  figure  7. 

Six  different  tower  sizes  were  used  in  the  analysis.  The 
effectiveness  model  underpredicts  the  tower  fan  energy  used  in  all 
cases.  The  largest  error,  13.7%,  is  found  at  tower  number  19  which 
denotes  a  Baltimore  Aircoil,  series  3000,  model  number  3269. 
Savings  in  computational  time  is  very  impressive  for  this  test 
variable.  Approximately  34  minutes  of  saved  time  is  seen  for  all 
cases,  roughly  an  80%  improvement  in  computational  run  time. 
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ACCURACY  OF  JOYCE  VS  EFFECTIVENESS 

FX)R  VARIOUS  TOWER  SIZES 
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figure  6.  Plot  of  Effectiveness  versus  Joyce  Data  for  Various 
Tower  Sizes. 

The  final  variable  manipulated  is  condenser  water  pump  size.  Pump 
sizes  between  750  GPM  and  1500  GPM  were  used  and  the  affect  on  both 
fan  and  pximp  energy  (KWH)  were  tabulated.  The  results  are 
promising,  producing  like  trends  for  both  the  effectiveness  model 
and  Joyce  method,  refer  to  figure  8.  Computational  run  time 
savings  are  also  favorable  for  this  parameter.  Figure  9  presents 
the  time  savings  for  the  new  program  as  compared  with  the  existing 
program  for  various  pump  sizes. 
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TIME  SAVINGS  AT  VARIOUS 

TOWER  SEES 
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Figur*  7.  Time  Savings  at  Various  Tower  Sizes 


Again  the  effectiveness  model  underpredicts  the  fan  and  pump  energy 
as  compared  to  the  Joyce  method  in  all  cases.  The  trends  for  both 
models,  however,  are  very  similar  giving  the  same  minimum  and 
maximum  point.  The  maximum  error  is  roughly  10%  and  occurs  at  a 
pump  size  of  1000  GPM.  Computer  run  time  savings  is  very 
impressive  at  approximately  34  minutes  or  81%  as  compared  to  the 
control  data  program. 

As  is  shown  in  all  the  figures,  the  trends  for  both  the 
effectiveness  model  and  the  Weber  model  in  Joyce's  program  are  very 
similar. 
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ACCURACY  OF  JOYCE  VS  EFFECTIVENESS 

PUMP  SEE 


Figure  8.  Plot  of  Effectiveness  versus  Joyce  Data  for  Various 
Pump  Sizes. 


In  order  to  improve  the  accuracy  of  the  effectiveness  method, 
crossflow  correction  factor  should  be  applied  and  the  water  loss 
from  the  cooling  tower  should  be  included.  The  affect  on  computer 
run  time,  however,  will  be  adverse  as  several  other  iteration  loops 
would  have  to  be  utilized.  Another  improvement  could  be  using  a 
more  accurate  curve  fit  for  the  saturation  specific  heat  C,.  This 
will  also  take  additional  computational  time  and  would  result  in 
only  a  minimal  improvement  in  accuracy. 
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The  results  of  this  project  are  very  encouraging.  For  a 
general  trend  and  a  first  tine  run  through  finding  the  optimun 
design  this  revised  effectiveness  model  fits  the  bill.  The 
reduction  in  computational  time  will  allow  an  engineer  flexibility 
in  investigating  various  systems  reaching  an  optimum  design  in  a 
shorter  time  frame. 
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EFFECTIVENESS  COMPUTER  PROGRAM 


17 


/ 


EFFECTIVENESS  PROGRAM  MODEL 


'  15  OCT  91 

'  THIS  PROGRAM  CALCULATES  THE  ANNUAL  ENERGY  COST,  AND  THE 

•  CAPITAL  CONSTRUCTION  COST  OF  THE  CONDENSER  SIDE  OF  A  CENTRIFUGAL 
'  CHILLER  PLANT.  VARIES  THE  CONDENSER  GPM  TO  MAINTAIN  OPT  GPM  FOR 
'  A  PARTICULAR  TEMP.  BIN  ASSUMING  1071  GPM  MAX.  USED  TO  COMPARE 

'  FIXED  SPEED  DESIGN  WITH  VARIABLE  SPEED. 

/ 

'  INPUT  PARAMETERS  t 

I 

•  1)  THE  CITY  (ATLANTA=1,  L.A.=2,  N.Y.=3) 

'  2)  THE  COOLING  TOWER  NO.  CORRESPONDING  TO  ONE  OF 

'  THE  COOLING  TOWERS  IN  THE  DOCUMENTATION. 

'  3)  THE  CONDENSER  WATER  PIPE  SIZE  IN  INCHES. 

'  4)  THE  COOLING  TOWER  CFM. 

'  5)  THE  NO.  OF  CONDENSER  PASSES  IN  THE  CHILLER. 

'  6)  THE  MAXIMUM  FLOWRATE  (GPM)  FOR  VARIABLE  SPEED 

'  PUMP. 

/ 

/  this  program  is  a  modification  of  THE  PROGRAM  ANUALFIX 

•  WRITTEN  BY  ERIC  WEBER  IN  CONJUNCTION  WITH  HIS  MASTERS 

'  THESIS  ENTITLED  "MODELING  AND  GENERALIZED  OPTIMIZATION  OF 

'  COMMERCIAL  BUILDING  CHILLER/ COOLING  TOWER  SYSTEMS” 

/ 

' _ UPDATES _ 

'  LAST  UPDATE  TO  ANUALFIX  BY  EDW  ON  23  OCT  1988 

'  MOD  TO  PREVENT  RECALC  OF  PUMP  HP  BY  CTJ  ON  3  SEPT  89 

'  MOD  TO  PREVENT  RECALC  OF  FAN  HP  BY  CTJ  ON  5  SEPT  89 

'  MOD  TO  TWR  CALC.  FOR  COOLING  WITH  FAN  OFF  BY  CTJ  ON  15  SEP  89 

•  MODI  TO  BISECTION  USE  OLD  CW  TEMP  FOR  BRACKETS  BY  CTJ  ON  29 

'  SEP  89 

'  MOD2  TO  BISECTION  TO  CHECK  SIGNS  ON  BRACKETS  BY  CTJ  ON  30 

'  SEP  89 

'  MOD3  TO  BISECTION  WARN  IF  SIGNS  ARE  SAME  BY  CTJ  ON  1  OCT  89 

'  MOD:  STD  WTHR  DATA  NOT  COMPR.  WTHR  DATA  BY  CTJ  ON  15  OCT  89 

"  MOD  TO  ADD  WTHR  DATA  FOR  N.Y.  AND  L.A.  BY  CTJ  ON  20  OCT  89 

'  MOD  TO  ADD  ENERGY  BILL  CALC.  FOR  3  CITIES  BY  CTJ  ON  23  OCT  89 
'  MOD  TO  ADD  CAPITAL  COST  FOR  3  CITIES  BY  CTJ  ON  25  NOCT  89 

'  MOD  TO  USE  BRENT  METHOD  VICE  OF  BISECTION  BY  CTJ  ON  27  OCT  89 
'  MOD  TO  ASSIGN  TCI  70  IF  PREVIOUS  HIGHER  TDB  ASSIGNED  TCI  TO 

'  70,  AND  CHECK  FOR  INCREASE  IN  DELTATCOND  BY  CTJ  ON  29  OCT  89 

'  MOD  TO  ADD  VARIABLE  SPEED  PUMP  BY  CTJ  ON  29  OCT  89 

'  UPDATE  TO  INCLUDE  THE  EFFECTIVENESS  MODEL  TO  CALCULATE  THE 

'  OUTPUT  TEMP  OF  THE  COOLING  TOWER  ON  15  OCT  91  BY  JDR 

f 

DECLARE  SUB  EFFECTIVE  (NTU! ,  TWIN!,  MDOTW! ,  MDOTA! ,  TWB! ,  TDB!, 
EXPTWO!,  TWOUT!) 


DECLARE  SUB  BRENT  (TWINl!,  TWIN2!,  DIFFll,  DIFF2 ! ,  TWOUTl ! , 
TWOUT2!,  NTU!,  fda!,  MDOTWl ,  MDOTA! ,  TDBi ,  TWB! ,  DELTATCONDl ! , 
CTR2%,  TCI!,  TC2!,  CITY!) 

DECLARE  SUB  SYSCOST  (PMPHP! ,  FANHP!,  TWRCOST! ,  diameter!,  CITY!, 
CAPCST ! ) 

DECLARE  SUB  ENRGYCST  (KWPK!(),  KWHTOT!(),  BILLDEM!(),  BILL!(), 
YRBILL!,  PMPKWH!(),  FANKWH!(),  CMPKWH!(),  CITY!) 

DECLARE  SUB  Hairm  (T2!,  Tl!,  W2!,  OMEGMIXTURE! ,  HAIRMIX!) 

DECLARE  SUB  Hairs  (TF! ,  WS! ,  HSAT!) 

DECLARE  SUB  FLOW  (GPM! ,  CFM! ,  TDB!,  MDOTW! ,  MDOTA!) 

DIM  HRSDAT(20),  KWPK(20) ,  KWHTOT(20) ,  BILLDEM(20) ,  BILL{20) , 
TWBR(20),  TDBR(20),  PMPKWH(20),  FANKWH(20),  CMPKWH(20) 

•  USER  INPUT 

'  -kit  It  It  ■kick  It  hit -kit  hit  it  hit  It  It  It  it  it  it  it  it  it  it 
/ 

INPUT  "INPUT  CITY  NUMBER  (1-3)  ",  CITY 

INPUT  "INPUT  TOWER  NUMBER  (1-19)  ",  i 

INPUT  "INPUT  DESIRED  CONDENSER  WATER  PIPE  SIZE  IN  INCHES  ", 
diameter 

INPUT  "INPUT  DESIRED  TOWER  CFM  ",  CFM 

INPUT  "INPUT  NUMBER  OF  CONDENSER  PASSES  ",  PASSES 

INPUT  "INPUT  THE  MAXIMUM  GPM  FOR  A  VARIABLE  SPEED  PUMP  ",  MAXGPM 

t  hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh 

"STORE  THE  START  TIME  FOR  THE  ITERATIONS 

I hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh 
/ 

STARTIME$  =  TIME$ 

f  hhhhhhhhhhhhhhhhhhhhhhh'ihhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh 

•  GET  THE  TOWER  CONSTANT  AND  THE  EXPONENT  FROM  DATA  FILE 

• hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh 
/ 


KWHSUM  =  0! 

TWR  =  i 

OPEN  "b: tower.dat"  FOR  INPUT  AS  #1 
DO  WHILE  i  >  0 

INPUT  #1,  TNUM,  EXPONENT,  CONSTANT,  GPML,  GPMH,  XXX,  FHPCA, 
FHPCB,  TWRCST 
i  =  i  -  1! 

LOOP 

CLOSE  #1 

9 

"  THIS  SECTION  SETS  UP  THE  DRY  BULB  TEMPERATURE  BINS 

>  hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh 
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FOR  k  =  1  TO  10 

TDBR(k)  =  102  -  ((k  -  1)  *  5) 

NEXT  k 
/ 

^  *****************  ***it*4t******************;k4t*4r4t**4;*it*4t*4tA******** 

'  THIS  SECTION  GETS  THE  WEATHER  DATA  FOR  THE  DRY  BULB  TEMP.  BINS  ' 
'  FOR  EACH  MONTH  —  JAN(i=l)  thru  Dec(i=12) 

/ 

FOR  i  =  1  TO  12 
/ 

/********************************** 

'  INITIALIZE  MONTHLY  KWH  AND  PEAKKW 

/********************************** 

KWHTOT(i)  =  0!  'MONTHLY  KWH 

KWPK(i)  =  01  'MAX  MONTHLY  KW  FOR  AT  LEAST  15%  OF  HRS 

'  IN  MONTH 

PMPKWH(i)  =  01 
FANKWH(i)  =  01 
CMPKWH(i)  =  01 

CWSTEMPCHK  =  0  'CHECK  VALUE  TO  INSURE  CWS  TEMP  DOESNT 

'  RISE  WITH  DECREASING  TDB,  AND  LOAD 
CTR%  =  0  'TRACKS  #  OF  TEMP  BINS  WITH  NON-0  TWB 

'  PER  MONTH 

HRSMNTH  -  01 
SELECT  CASE  CITY 
CASE  1 

OPEN  "b; atlwth.dat”  FOR  INPUT  AS  #3 
CASE  2 

OPEN  "b;  lawth.dat*'  FOR  INPUT  AS  jf3 
CASE  3 

OPEN  ''b:nywth.dat''  FOR  INPUT  AS  #3 
END  SELECT 
k  =  i 

DO  WHILE  k  >  01 
FOR  m  =  1  TO  10 

INPUT  #3,  HRSDAT(m),  HOURSSUM,  TWBR(m) 

NEXT  m 
k  =  k  -  1 
LOOP 
CLOSE  i3 

/ 

'****************************** 

'SELECT  THE  WEATHER  DATA  BIN 
/****************************** 
t 

FOR  j  =  1  TO  10 
TWB  =  TWBR(j) 

EXPTWO  =  TWB 


'WET  BULB  TEMP 

'EXPECTED  TEMP  WATER  FROM  TOWER 
'  FOR  EFF  MODEL 
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TDB  =  TDBR(j) 
HOURS  =  HRSDAT(j) 


'DRY  BULB  TEMP 

'MONTHLY  HOURS  AT  TDB  &  TWB  FROM 
'  6AM  TO  6PM 

IF  TWB  =  0  THEN  GOTO  GETNEXTWEATHER 
CTR%  =  CTR%  +  1 


CTR1% 

CTR2% 

CTR3% 

CTR4% 


=  0 


=  0 


=  0 


=  0 


'TRACKS  #  OF  ITTERATIONS  THRU 
'  RESTARTA 

'TRACKS  #  OF  ITTERATIONS  THRU 
'  BRENT  MTHD 

ITTERATIONS  THRU 


'TRACKS  #  OF 
'  TRYAGAIN 
'TRACKS  #  OF 
'  CYCLEFAN 


ITTERATIONS  THRU 


> ************************************** 
'CALCULATION  OF  THE  EVAPORATING  TEMP. 
'*********************************11**** 


EVAPGPM  =  1200 
TE2  =  44 


> ****************************** 

•  EVAP  LOAD  ALGORITHM 

I  ****************************** 

/ 

GPM  ==  -1242  +  (93.90529  *  TDB)  -  (1.45875  *  TDB  ^  2)  +  (.007462 

•  TDB  ''  3) 

IF  GPM  <=  775  THEN  GPM  =  750 
tons  =  -531.578  +  10.526  *  TDB 
IF  tons  >500  THEN 
tons  =  500 
END  IF 

IF  tons  <  100  THEN 
tons  =  100 
END  IF 

QDOTEVAP  =  tons  *  12000 

DELTATEVAP  =  QDOTEVAP  /  (EVAPGPM  *  499.27) 

UAEVAP  =  1000  *  (134.11603#  +  1284.3418#  *  (tons  /  500)  - 

640.54974#  *  ((tons  /  500)  ^  2)) 

LMTDevap  =  QDOTEVAP  /  UAEVAP 

TE  =  (-TE2  -  DELTATEVAP  +  TE2  *  EXP ( DELTATEVAP  /  LMTDevap))  / 
(EXP (DELTATEVAP  /  LMTDevap)  -  l!) 

•  *********************************1,1, 

•  ESTIMATE  OF  THE  KW  FOR  COMPRESSOR 
************************************ 

/ 

KW  =  30  +  .61  *  (tons) 


'FLOWRATE  THRU  EVAPORATOR 
'FIXED  CHILLED  WATER  SUPPLY 

/  TPMT5. 
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TWINl  =50.5 
TWIN2  =  110 

t 

r***************************************************** 

'CALCULATE  THE  HEAT  BURDEN  FOR  THE  CHILLER  CONDENSER 

f********ic********1i**1i'k1i1t1i1t1fk1i****1c*1fk****ii**1t1t1c*1t***1t 


RESTARTA: 

QDOTCOND  =  QDOTEVAP  +  (KW  *  3413) 

DELTATCONDl  =  QDOTCOND  /  (GPM  *  499.27) 

CALL  FLOW (GPM,  CFM,  TDB,  MDOTW,  MDOTA) 

fda  =  MDOTW  *  CONSTANT  *  ((MDOTW  /  MDOTA)  ''  EXPONENT) 

NTU  =  CONSTANT  *  (  (MDOTW  /  MDOTA)  (EXPONENT  +  1)  ) 

/ 

'  BRACKETING  FOR  BRENT  METHOD 

f******************************** 

/ 

CALL  EFFECTIVE (NTU,  TWINl,  MDOTW,  MDOTA,  TWB,  TDB,  EXPTWO, 

TWOUTl) 

DELTAl  =  TWINl  -  TWOUTl 
DIFFl  =  DELTATCONDl  -  DELTAl 

CALL  EFFECTIVE (NTU,  TWIN2 ,  MDOTW,  MDOTA,  TWB,  TDB,  EXPTWO, 

TWOUT2 ) 

DELTA2  =  TWIN2  -  TWOUT2 
DIFF2  =  DELTATCONDl  -  DELTA2 
IF  SGN (DIFFl)  =  SGN(DIFF2)  THEN 
TWINl  =  TWINl  -  5 
TWIN2  =  TWIN2  +  5 
SELECT  CASE  CITY 
CASE  1 

OPEN  ”B:ATLEFF.DAT''  FOR  APPEND  AS  #4 
CASE  2 

OPEN  "BzLAEFF.DAT"  FOR  APPEND  AS  #4 
CASE  3 

OPEN  "B.-NYEFF.DAT”  FOR  APPEND  AS  #4 
END  SELECT 

PRINT  #4,  "ROOT  NOT  BRACKETED  FOR  BRENT.  PROG  ASSIGNING  NEW 
BRACKETS" 

PRINT  #4,  "MONTH=  ";  i,  "BIN  =",'  TDB 
CLOSE  #4 

PRINT  "ROOT  NOT  BRACKETED  FOR  BRENT.  PROG  ASSIGNING  NEW 
BRACKETS" 

PRINT  "MONTH=  "  ,*  i,  "BIN  =";  TDB 
GOTO  RESTARTA: 

END  IF 

CTR1%  =  CTR1%  +  1 

CALL  BRENT (TWINl,  TWIN2 ,  DIFFl,  DIFF2,  TWOUTl,  TWOUT2 ,  NTU,  fda, 
MDOTW,  MDOTA,  TDB,  TWB,  DELTATCONDl,  CTR2%,  TCI,  TC2,  CITY) 

DELTA3  =  TWIN2  -  TWOUT2 
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/ 


'  DETERMINATION  OF  CONDENSER  PERFORMANCE  (UA,Tc) 
'************************************************* 


GPMCOND  =  GPM 

UAcond  =  608.5483400000001#  +  130.64856#  *  (tons  /  500)  - 
192.20566#  *  ((tons  /  500)  ^  2) 

UAcond  =  UAcond  +  1203.80505#  *  (GPM  /  (PASSES  *  1500))  - 
557.4549#  *  ((GPM  /  (PASSES  *  1500))  "  2) 

UAcond  =  UAcond  *  1000 
LMTDCOND  =  QDOTCOND  /  UAcond 

TC  =  (TCI  -  (TCI  +  (TC2  -  TCI))  *  EXP((TC2  -  TCI)  /  LMTDCOND))  / 

(1!  -  EXP((TC2  -  TCI)  /  LMTDCOND)) 

/ 

'  *********************************■)(***************************** 

'DETERMINE  MODEL  KW  USING  THE  CARNOT  KW/TON  AND  TWO  EFFICIENCIES 
^ *************************************************************** 
t 

KWCARNOT  =  (((TC  +  460)  /  (TE  +  460))  -  l!)  *  3.517 
NfTON  =  4,5869  *  (tons  /  500)  -  8.16536  *  ((tons  /  500)  "  2)  + 
6.65014  *  ((tons  /  500)  "  3)  -  2.07167  *  ((tons  /  500)  "  4) 

NfISENT  =  .516061  +  ,31889  *  ((KWCARNOT  -  .3)  /  .2)  -  .14911#  * 
(((KWCARNOT  -  .3)  /  ,2)  "  2) 

KWTONMODEL  =  KWCARNOT  /  (NfISENT  *  NfTON) 

KWTONGUESS  =  KW  /  tons 

/ 

^ ********************************************** 

'  COMPARE  MODEL  KW/TON  TO  THE  GUESSED  KW/TON 
^********************************************** 

IF  ABS (KWTONMODEL  -  KWTONGUESS)  >  .0001  THEN 
KW  =  KWTONMODEL  *  tons 
GOTO  RESTARTA 
END  IF 
TCIA  =  TCI 
TC2A  =  TC2 

/ 

^ **************************************************************** 
'  IF  THE  ENTERING  COND.  WATER  TEMP.  IS  BELOW  80  THEN  FIX  IT  AT  80 
'  AND  CYCLE  THE  COOLING  TOWER  FAN. 
'**************************************************************** 
/ 

IF  TCI  <  80  THEN 
GOTO  CYCLEFAN 
END  IF 

IF  CTR%  >  1  THEN 
GOTO  CALCTOTALKW 
END  IF 

FANHP  =  EXP(FHPCA  +  ((LOG(CFM)  -  10.5)  /  21)  *  FHPCB) 

FanKW  =  ((FANHP  *  .7457)  /  .85) 

GOTO  CALCTOTALKW 
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•  THIS  SECTION  WILL  HOLD  CONDENSER  WATER  INPUT  TEMPERATURE  TO 
'  80  F  AND  CYCLE  THE  TOWER  FAN  USING  THE  TOWER  TEMPERATURE  DROP 

'  VERSUS  THE  CONDENSER  TEMPERATURE  DROP 

9 

^  'k'k •kit it ’k'k'k’k 
9 

CYCLEFAN: 

TRYAGAIN; 

CTR3%  =  CTR3%  +  1 

QDOTCOND  =  QDOTEVAP  +  (KW  *  3413) 

DELTATCOND  *=  QDOTCOND  /  (GPM  *  499.27) 

CALL  FLOW (GPM,  CFM,  TDB,  MDOTW,  MDOTA) 

fda  =  MDOTW  *  CONSTANT  *  ( (MDOTW  /  MDOTA)  ^  EXPONENT) 

NTU  =  CONSTANT  *  ((MDOTW  /  MDOTA)  ^  (EXPONENT  +1)) 

TCI  =  80 

TC2  =  TCI  +  DELTATCOND 

CALL  EFFECTIVE (NTU,  TC2 ,  MDOTW,  MDOTA,  TWB,  TDB,  EXPTWO,  TWOUTl) 
DELTATTOWER  =  TC2  -  TWOUTl 

UAcond  =  608.5483400000001#  +  130.64856#  *  (tons  /  500)  - 

192.20566#  *  ((tons  /  500)  ^  2) 

UAcond  =  UAcond  +  1203.80505#  *  (GPM  /  (PASSES  *  1500))  - 

557.4549#  *  ((GPM  /  (PASSES  *  1500))  ^  2) 

UAcond  =  UAcond  *  1000 
LMTDCOND  =  QDOTCOND  /  UACond 

TC  =  (TCI  -  (TCI  +  (TC2  -  TCI))  *  EXP((TC2  -  TCI)  /  LMTDCOND))  / 
(1!  -  EXP((TC2  -  TCI)  /  LMTDCOND)) 

KWCARNOT  =  (((TC  +  460)  /  (TE  +  460))  -  l!)  *  3.517 

NfTON  =  4.5869  *  (tons  /  500)  -  8.16536  *  ((tons  /  500)  "  2)  + 

6.65014  *  ((tons  /  500)  "  3)  -  2.07167  *  ((tons  /  500)  "  4) 

NfISENT  =  .516061  +  .31889  *  ((KWCARNOT  -  .3)  /  .2)  -  .14911#  * 
(((KWCARNOT  -  .3)  /  .2)  ^  2) 

KWTONMODEL  =  KWCARNOT  /  (NfISENT  *  NfTON) 

KWTONGUESS  =  KW  /  tons 

IF  ABS (KWTONMODEL  -  KWTONGUESS)  >  .0001  THEN 
KW  =  KWTONMODEL  *  tons 
GOTO  TRYAGAIN 
END  IF 

FANHP  =  EXP(FHPCA  +  ((LOG(CFM)  -  10.5)  /  2!)  *  FHPCB) 

TOWERDUTY  =  (DELTATCOND  -  (.1  *  DELTATTOWER))  /  (DELTATTOWER  * 
.9) 

FanKW  =  ((FANHP  *  .7457)  /  .85)  *  TOWERDUTY 
CTR4%  =  CTR4%  +  1 

t 

>**************************************************************** 

'CALCULATION  OF  THE  POWER  REQUIRED  BY  THE  CONDENSER  WATER  PUMP 
f****************************** ********************************** 
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CALCTOTALKW: 

PipeD  =  diameter  /  12 

Condenservelocity  =  (GPM  *  .002228)  *  (PASSES  /  1.15241) 
Pipevelocity  «  (GPM  *  .002228)  /  ((3.1416  *  (PipeD  ^  2))  /  4) 
REd  =  (Pipevelocity  *  PipeD)  /  8.64E-06 
IF  REd  <  2300  THEN 
pipefrict  =  64  /  REd 
ELSE 

'  FROM  EQU.  6.64a  IN  FLUID  MECHANICS  BY  WHITE 

'  It  It  It  h  It  It  It  It  it  It  It  It  it  It  It  h  it  it  it  It  It  It  h  •kit  it  it  it  It  It  It  it  It  it  It  hit  It  It  It  It  It 

AA  =  (6.9  /  REd)  +  ((.00015  /  (PipeD  *  3.7))  "  1.11) 

BB  =  LOG(AA)  /  LOG(lOl) 
pipefrict  =  (1!  /  (-1.8  *  BB) )  2 

END  IF 

PHtower  =11 
SELECT  CASE  PASSES 
CASE  1 

PHcond  =  EXP(-13. 95945  +  20.21629  *  (LOG(GPM)  /  10)) 

CASE  2 

PHcond  =  EXP(-10. 96759  +  19.02989  *  (LOG(GPM)  /  10)) 

CASE  3 

PHcond  =  EXP(-8. 96627  +  17.76877  *  (LOG(GPM)  /  10)) 

CASE  4 

PHcond  =  EXP (-8. 76329  +  18.62906  *  (LOG (GPM)  /  10!)) 

END  SELECT 

PipeL  =  200 

PipefL  =  450  *  PipeD 

PipestrainL  =  250  *  PipeD 

Ltotal  =  PipeL  +  PipefL  +  PipestrainL 

PHpipe  =  ((Pipevelocity  ''  2)  /  64.4)  *  ((pipefrict  *  Ltotal)  / 
PipeD) 

PHTOTAL  =  PHcond  +  PHtower  +  PHpipe 
PERGPM  =  (GPM  /  MAXGPM)  *  100 

PMPEFF  =  10.57993  +  (2.436639  *  PERGPM)  -  (.01966  *  (PERGPM  ''2)) 
VFDEFF  =  4.983855  +  (1.846339  *  PERGPM)  -  (.01029  *  (PERGPM  ^  2)) 
PMPHP  =  (GPM  *  PHTOTAL)  /  (3960  *  PMPEFF  /  100) 

PumpKW  =  (PMPHP  *  .7457)  /  (VFDEFF  /  100) 

t 

f  hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh 

•  CALCULATION  OF  THE  TOTAL  POWER  REQUIRED  BY  THE  SYSTEM 
'hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh 


TOTALKW  =  PumpKW  +  FanKW  +  KW  'TOTAL  SYSTEM  KW 
KWHBIN  =  TOTALKW  *  HOURS  'ENERGY  CONSUMPTION  OF  BIN 

KWHTOT(i)  =  KWHTOT(i)  +  KWHBIN  'ENERGY  CONSUMPTION  OF  MONTH 
KWHSUM  =  KWHSUM  +  KWHBIN  'ENERGY  CONSUMPTION  OVER  YEAR 

PKHOURS  =  .15  *  HOURSSUM 

IF  (HOURS  >=  PKHOURS)  AND  (TOTALKW  >  KWPK(i))  THEN  KWPK(i)  = 
TOTALKW 
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HRSMNTH  *  HRSMNTH  +  HOURS 

PMPKWH(i)  =  PMPKWH(i)  +  (PumpKW  *  HOURS) 

FANKWH(i)  =  FANKWH(i)  +  (FanKW  *  HOURS) 

CMPKWH(i)  =  CMPKWH(i)  +  (KW  *  HOURS) 

f 

^  it  It  it  it  icit  it  it  it  it  It  ic  it  It  it  it  it  It  it  it  It  it  "kit  It  It  it  it  if -kit  It  it  4t  if 

"STORE  THE  END  TIME  FOR  ITERATIONS 

'A********************************** 

/ 

ENDTIME$  =  TIME$ 

/ 

"  PRINTING  OF  THE  RESULTS 
/ 

SELECT  CASE  CITY 
CASE  1 

OPEN  "B: ATLEFF.DAT”  FOR  APPEND  AS  #4 
CASE  2 

OPEN  ''B:LAEFF.DAT''  FOR  APPEND  AS  #4 
CASE  3 

OPEN  "BzNYEFF.DAT''  FOR  APPEND  AS  /4 
END  SELECT 

PRINT  #4,  "TIME  START  =  ";  STARTIME$;  "TIME  END  =  ",  ENDTIME$ 
PRINT  #4,  "Ctrl*”;  CTRl%,  "Ctr2  =  CTR2%,  ••  Ctr3  =  CTR3%; 
PRINT  #4,  ''I*”;  i,  ”J=";  j 

PRINT  #4,  ''TDB(DEG  F)  TWB(DEG  F)  HOURS  TOWER  *• 

PRINT  #4,  ;  TDB,  TWB,  HOURS,  TWR 

PRINT  #4,  "GPM  %GPM  CFM  TONS" 

PRINT  #4,  ;  GPM,  PERGPM,  CFM,  tons 
IF  TCI  >80  THEN 

PRINT  #4,  ''CWS(DEG  F)  CWR(DEG  F)  DELTATCOND  DELTATTOWER" 

PRINT  #4,  ;  TCIA,  TC2A,  DELTATCONDl,  DELTA3 

ELSE 

PRINT  #4,  ••CWS(DEG  F)  CWR(DEG  F)  DELTATCOND  DELTATTOWER 

TOWERDUTY" 

PRINT  #4,  ;  TCI,  TC2 ,  DELTATCOND,  DELTATTOWER,  TOWERDUTY 
END  IF 

PRINT  #4,  “PIPE  D(FT)  Pipevelocity (FT/SEC)  REd  pipe 

frict  ” 

PRINT  #4,  ;  PipeD,  Pipevelocity,  REd,  pipefrict 
PRINT  #4,  "pipeL  total  PHpipe  PHcond  PHtower 

PHtotal(FT  Water)” 

PRINT  #4,  ;  Ltotal,  PHpipe,  PHcond,  PHtower,  PHTOTAL 
PRINT  #4,  ”PASSES  Condenservelocity” 

PRINT  #4,  ;  PASSES,  Condenservelocity 

9 

PRINT  #4,  ”UACOND  LMTDCOND  NfISENT  NFTON 

COND.TEMP" 

PRINT  #4,  ;  UAcond,  LMTDCOND,  NfISENT,  NfTON,  TC 

PRINT  #4,  "KWTONGUESS  KWTONMODEL  KWHBIN  KWHTOT" 
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PRINT 

#4, 

;  KWTONGUESS, 

KWTONMODEL,  KWHBIN,  KWHTOT(i) 

PRINT 

#4, 

"HOURS 

HOURSSUM" 

PRINT 

#4, 

;  HOURS,  HOURSSUM 

PRINT 

#4, 

"CHILLERKW 

PUMPKW  FANKW  TOTALKW 

PRINT 

#4, 

;  KW,  PumpKW, 

FanKW,  TOTALKW 

PRINT 

#4, 

II II 

PRINT 

#4, 

It  11 

PRINT 

#4, 

nn 

PRINT 

#4, 

II II 

CLOSE 

#4 

OVER: 

GETNEXTWEATHER : 

f  ******************************* 

'GET  THE  NEXT  WEATHER  DATA  BIN 

f  ******************************* 

NEXT  j 

/******************************* 

'START  CALC  FOR  THE  NEXT  MONTH 

t ******************************** 

NEXT  i 
/ 

CALL  ENRGYCST(KWPK() ,  KWHTOT() ,  BILLDEM() ,  BILL(),  YRBILL, 
PMPKWH ( ) ,  FANKWH ( ) ,  CMPKWH ( ) ,  CITY) 

CALL  SYSCOST(PMPHP,  FANHP,  TWRCST,  diameter,  CITY,  CAPCST) 

/ 

END 


SUB  BRENT  (TWINl,  TWIN2 ,  DIFFl,  DIFF2,  TWOUTl,  TWOUT2 ,  NTU,  fda, 
MDOTW,  MDOTA,  TDB,  TWB,  DELTATCONDl,  CTR2%,  TCI,  TC2,  CITY) 

ITMAX  =  50 
TOL  =  .0001 
EPS  =  1  *  10  "  -7 
DIFF3  =  DIFF2 
TWOUT3  =  TWOUT2 
FOR  ITER  =  1  TO  ITMAX 
CTR2%  =  CTR2%  +  1 
IF  (DIFF2  *  DIFF3)  >  0  THEN 
TWIN3  =  TWINl 
DIFF3  =  DIFFl 
TWOUT3  =  TWOUTl 
TWIN4  =  TWIN2  -  TWINl 
TWINS  =  TWIN4 
END  IF 

IF  ABS(DIFF3)  <  ABS(DIFF2)  THEN 
TWINl  =  TWIN2 
TWIN2  =  TWIN3 
TWIN3  =  TWINl 
DIFFl  =  DIFF2 
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TWOUTl  =  TWOUT2 
DIFF2  =  DIFF3 
TWOUT2  =  TWOUT3 
DIFF3  =  DIFFl 
TWOUT3  =  TWOUTl 
END  IF 

TOLl  =  (2  *  EPS  *  ABS(TWIN2))  +  (.5  *  TOL) 

XM  *  .5  *  (TWIN3  -  TWIN2) 

IF  (ABS(XM)  <=  TOLl)  OR  (DIFF2  =  0)  THEN 
TC2  =  TWIN2 
TCI  =  TWOUT2 
EXIT  SUB 
END  IF 

IF  (ABS(TWIN5)  >=  TOLl)  AND  (ABS(DIFFl)  >  ABS(DIFF2))  THEN 
S  =  DIFF2  /  DIFFl 
IF  TWINl  =  TWIN3  THEN 
P  =  2  *  XM  *  S 
Q  =  1  -  S 
ELSE 

Q  =  DIFFl  /  DIFF3 
R  =  DIFF2  /  DIFF3 

P  =  S  *  (2  *  XM  *  Q  *  (Q  -  R)  -  (B  -  A)  *  (R  -  1)  ) 

Q  =  (Q  -  1)  *  (R  -  1)  *  (S  -  1) 

END  IF 

IF  P  >  0  THEN  Q  *  -Q 
P  =  ABS(P) 

MINI  =  3  *  XM  *  Q  -  ABS(T0L1  *  Q) 

MIN2  =  ABS (TWINS  *  Q) 

IF  MINI  <  MIN2  THEN 
MIN  =  MINI 
ELSE 

MIN  =  MIN2 
END  IF 

IF  (2  *  P)  <  MIN  THEN 
TWINS  =  TWIN4 
TWIN4  =  P  /  Q 
ELSE 

TWIN4  =  XM 
TWINS  «  TWIN4 
END  IF 
ELSE 

TWIN4  =  XM 
TWINS  =  TWIN4 
END  IF 

TWINl  =  TWIN2 
DIFFl  =  DIFF2 
TWOUTl  =  TWOUT2 
IF  (ABS(TWIN4)  >  TOLl)  THEN 
TWIN2  =  TWIN2  +  TWIN4 
ELSE 

SELECT  CASE  XM 
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CASE  IS  >  0 
SIGN  =  1 
CASE  0 
SIGN  =  1 
CASE  IS  <  0 
SIGN  =  -1 
END  SELECT 
T0L2  =  TOLl  *  SIGN 
TWIN2  *  TWIN2  +  TOL2 
END  IF 

CALL  EFFECTIVE (NTU,  TWIN2,  MDOTW,  MDOTA,  TWB,  TDB,  TWB,  TWOUT2) 
DELTA2  =  TWIN2  -  TWOUT2 
DIFF2  =  DELTATCONDl  -  DELTA2 
NEXT  ITER 

SELECT  CASE  CITY 
CASE  1 

OPEN  •*B:ATLEFF.DAT''  FOR  APPEND  AS  #1 
CASE  2 

OPEN  ''BiLAEFF.DAT''  FOR  APPEND  AS  #1 
CASE  3 

OPEN  ''BrNYEFF.DAT”  FOR  APPEND  AS  #1 
END  SELECT 

PRINT  #1,  "BRENT  METHOD  EXCEEDED  MAX  ITERATIONS" 

TC2  =  TWIN2 
TCI  =  TWOUT2 
CLOSE  #1 

PRINT  "BRENT  METHOD  EXCEEDED  MAX  ITERATIONS" 

PRINT  "MONTH  =" ;  i,  "BIN  =";  TDB 
END  SUB 


SUB  EFFECTIVE  (NTU,  TWIN,  MDOTW,  MDOTA,  TWB,  TDB,  EXPTWO,  TW2) 
DELTAS  =10 
Ttemp  =  EXPTWO 

DO  WHILE  ABS (DELTAS)  >  .0001 

/ 

'Calculate  the  inlet  and  outlet  saturation  enthalpy  of  the 

'saturated  air. 

/ 

CALL  Hairs (Ttemp,  WSATl,  hswo) 

CALL  Hairs (TWIN,  WSATl,  hswi) 

/ 

'Calculate  the  saturation  specific  heat  and  the  m*  parameter  as 

'defined  by  Braun. 

/ 

Cs  =  (hswi  -  hswo)  /  (TWIN  -  Ttemp) 
m  =  MDOTA  /  (MDOTW  *  (1  /  Cs) ) 

eff  =  (1  -  EXP(-NTU  *  (1  -  m)))  /  (1  -  m  *  EXP(-NTU  *  (1  -  m))) 
'Calculate  the  humidity  ratio  of  the  air-vapor  mixture. 
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CALL  Hairs (TWB,  WSAT,  H) 

CALL  Hainn(TWB,  TDB,  WSAT,  OMEGl,  hai) 
hao  =  hai  +  eff  *  (hswi  -  hai) 

Tref  *  32 

TW2  =  ( (MDOTW  *  (TWIN  -  Tref)  -  MDOTA  *  (hao  -  hai) )  /  (MDOTW  *  1) ) 
+  Tref 

DELTA3  =  TW2  -  Ttemp 
Ttemp  =  TW2 
LOOP 
END  SUB 


SUB  ENRGYCST  (KWPK() ,  KWHTOT(),  BILLDEM(),  BILL() ,  YRBILL, 
PMPKWHO,  FANKWHO,  CMPKWH(),  CITY) 

TOTKWH  =  0 

FOR  1  =  1  TO  12 

TOTKWH  =  TOTKWH  +  KWHTOT(l) 

KWHTOT(l)  =  KWHTOT(l)  +  360585 
KWPK(l)  =  KWPK(l)  +  946 
NEXT  1 

SELECT  CASE  CITY 
CASE  1 

GOTO  ATLANTA 
CASE  2 

GOTO  LOSANGELES 
CASE  3 

GOTO  NEWYORK 
END  SELECT 
ATLANTA: 

f1t1i***1i*1fk*1t*1i1i*1t*1i***1i1i1fk'k**1flfii*1i**1Ht*1fk**1t***1t**'it*****ifk**it* 

•  CALCULATION  OF  MONTHLY  ENERGY  COST  IN  ATLANTA 

•  Georgia  Power  schedule  PL  8_GS 

•  General  Commercial  Service 

'  Current  as  of  June  1989 


KWDEM  =  5 
FOR  i  =  6  TO  9 

KWSMDEM  =  .95  *  KWPK(i) 

IF  KWSMDEM  >  KWDEM  THEN  KWDEM  =  KWSMDEM 
NEXT  i 

FOR  i  =  1  TO  12 

IF  (i  >  5)  AND  (i  <  10)  THEN  GOTO  ATLNOTWNTR 
KWWNTR  =  .6  *  KWPK(i) 

IF  KWWNTR  >  KWDEM  THEN  KWDEM  =  KWWNTR 
ATLNOTWNTR; 

NEXT  i 

FOR  i  =  6  TO  9 

IF  KWPK(i)  >  KWDEM  THEN 
BILLDEM(i)  =  KWPK(i) 

ELSE 
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BILLDEM(i)  =  KWDEM 
END  IF 
NEXT  i 

FOR  i  =  1  TO  5 

BILLDEM(i)  =  KWDEM 
NEXT  i 

FOR  i  =  10  TO  12 

BILLDEM(i)  =  KWDEM 
NEXT  i 

f*************************** **************************** 

•  CALCULATION  OF  THE  MONTHLY  ENERGY  BILL  FOR  ATLANTA 

/******************************************************* 


YRBILL  =  0 
PMPENGY  =  0 
FANENGY  =  0 
CMPENGY  =  0 
FOR  i  =  1  TO  12 

BLOKHRS  =  KWHTOT(i)  /  BILLDEM(i) 
IF  (BLOKHRS  <=  200)  THEN 
GOTO  BASEI 

ELSEIF  (BLOKHRS  100)  THEN 
GOTO  BASEI I 

ELSEIF  (BLOKH^<r  <=  600)  THEN 

GOTO  BASEI IT 

ELSE 

GOTO  BAS El V 
END  IF 

/ 


BASEI : 

IF  KWHTOT(i)  <=  3000  THEN 

COSTl  =  (.10185  *  KWHTOT(i))  + 
ELSEIF  KWHTOT(i)  <=  10000  THEN 
COSTl  =  (.09282  *  KWHTOT(i))  + 
ELSEIF  KWHTOT(i)  <=  200000  THEN 
COSTl  =  (.07918  *  KWHTOT(i))  + 
ELSE 

COSTl  =  (.05972  *  KWHTOT(i))  + 
END  IF 

GOTO  ATLSUMBILL 


13.5 

40.59 

176.99 

4068.99 


BASEII: 

IF  (BILLDEM(i)  >=  5)  AND  (BILLDEM(i)  <=  15)  THEN 

COSTl  =  (18.314  *  BILLDEM(i))  +  (.01028  *  KWHTOT(i))  +13.5 
ELSEIF  (BILLDEM(i)  <=  50)  THEN 

COSTl  =  (16.508  *  BILLDEM(i))  +  (.01028  *  KWHTOT(i))  +  40.59 
ELSEIF  (BILLDEM(i)  <=  1000)  THEN 

COSTl  =  (13.78  *  BILLDEM(i))  +  (.01028  *  KWHTOT(i))  +  176.99 
ELSE 

COSTl  =  (9.888  *  BILLDEM(i))  +  (.01028  *  KWHTOT(i))  +  4068.99 
END  IF 

GOTO  ATLSUMBILL 
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BASEIIX  z 

IF  (BILLDEM(i)  >=  5)  AND  (BILLDEM(i)  <=  15)  THEN 

COSTl  =  (19.938  *  BILLDEM(i))  +  (.00622  *  KWHTOT(i))  +13.5 
ELSEIF  BILLDEM(i)  <«  50  THEN 

COSTl  =  (18.132  *  BILLDEM(i))  +  (.00622  *  KWHTOT(i))  +  40.59 
ELSEIF  BILLDEM(i)  <=  1000  THEN 

COSTl  =  (15.404  *  BILLDEM(i))  +  (.00622  *  KWHTOT(i))  +  176.99 
ELSE 

COSTl  =  (11.512  *  BILLDEM(i))  +  (.00622  *  KWHTOT(i))  +  4068.99 
END  IF 

GOTO  ATLSUMBILL 

f 

BASEIV; 

IF  (BILLDEM(i)  >=  5)  AND  (BILLDEM(i)  <=  15)  THEN 

COSTl  =  (20.46  *  BILLDEM(i))  +  (.00535  *  KWHTOT(i))  +13.5 
ELSEIF  BILLDEM(i)  <=50  THEN 

COSTl  =  (18.654  *  BILLDEM(i))  +  (.00535  *  KWHTOT(i))  +  40.59 
ELSEIF  BILLDEM(i)  <=  1000  THEN 

COSTl  =  (15.926  *  BILLDEM(i))  +  (.00535  *  KWHTOT(i))  +  176.99 
ELSE 

COSTl  =  (12.034  *  BILLDEM(i))  +  (.00535  *  KWHTOT(i))  +  4068.99 
END  IF 

*******  CALC.  OF  THE  TOTAL  BILL  ******** 

ATLSUMBILL; 

FCADJ  =  .016896  *  KWHTOT(i) 

TAX  =  .05  *  (COSTl  +  FCADJ) 

BILL(i)  =  COSTl  +  FCADJ  +  TAX 
^*****CHECK  FOR  MIN.  BILL************ 

CHKKW  =  BILLDEM(i)  -  30 
IF  CHKKW  >  0  THEN 

MINBILL  =  (13.5  +  (7.5  *  CHKKW)  +  FCADJ)  *  1.05 
ELSE 

MINBILL  =  (13.5  +  FCADJ)  *  1.05 
END  IF 

IF  MINBILL  >=  BILL(i)  THEN  BILL(i)  =  MINBILL 
YRBILL  =  YRBILL  +  BILL(i) 

PMPENGY  =  PMPENGY  +  PMPKWH(i) 

FANENGY  =  FANENGY  +  FANKWH(i) 

CMPENGY  =  CMPENGY  +  CMPKWH(i) 

NEXT  i 


FOR  j  =  1  TO  12 

OPEN  "B; ATLEFF.DAT"  FOR  APPEND  AS 
PRINT  #4,  "THE  MAX  MONTHLY  KW  =", 
"THE  BILLING  DEMAND  =", 
"THE  MONTHLY  KWH  =" , 
"THE  MONTHLY  BILL  =  " 


PRINT  #4, 
PRINT  #4, 
PRINT  #4, 


PRINT  #4,  "" 
PRINT  #4,  "" 


#4 

KWPK(j) 

BILLDEM(j) 

KWHTOT(j) 

,  BILL(j) 


IF  j  =  12  THEN 
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PRINT  #4, 
PRINT  #4, 
PRINT  #4, 
PRINT  #4, 
PRINT  #4, 
PRINT  #4, 
PRINT  #4, 
END  IF 
CLOSE  #4 
NEXT  j 
EXIT  SUB 


••THE  YEARLY  TOTAL  ($)  =•*,  YRBILL 

••THE  PUMP  ENERGY  (KWH)  =•• ,  PMPENGY 
••THE  FAN  ENERGY  (KWH)  *••,  FANENGY 
••THE  CHILLER  ENERGY  (KWH)  =••,  CMPENGY 

••THE  TOTALSYSTEM  ENERGY  (KWH)  =•• ,  TOTKWH 
n  11 

n  II 


LOSANGELES : 

'**************************************************************** 
'  CALCULATION  OF  MONTHLY  ENERGY  COSTS  IN  LOS  ANGELES 

'  Souther  California  Edison  Schedule  GS-2 

'  General  Commercial  -  Large 

'  Eff.  Date  1  June  1988 

'  Current  as  of  June  1989 


'  **  IN  LA  THE  MONTHLY  DEMAND  IS  STRICTLY  DET.  BY  THE  CURRENT 

•  MONTH'S  PERFORMANCE.  THEREFORE  THE  DEMAND  COST  IS  AS  FOLLOWS: 

I  •k'k'k-k'kifk'klt'itlfklfklfk'klfklfklfkltlfkifklfklfklfklfkltlt'klfk'kle'kltlfk'k'klfk'klfk'k'k'kitlt'kltlfklt 

YRBILL  *  0 
PMPENGY  =  0 
FANENGY  =  0 
CMPENGY  =  0 
FOR  i  =  1  TO  12 

BILLDEM(i)  =  KWPK(i) 

DEMCOST  =  0 

IF  (i  >  5)  AND  (i  <  10)  GOTO  LASUMMER 
MIN  =  31.1 

DEMCOST  =  2.7  *  BILLDEM(i) 

IF  KWHTOT(i)  <=  180000  THEN 

ENGCOST  =  .07197  *  KWHTOT(i) 

ELSEIF  KWHTOT(i)  >  180000  THEN 

ENGCOST  =  12954.6  +  ((KWHTOT(i)  -  180000)  *  .04973) 

ELSE 
END  IF 

GOTO  LASUMBILL 
LASUMMER: 

MIN  =  31.1 

DEMCOST  =  8.6  *  BILLDEM(i) 

IF  KWHTOT(i)  <=  180000  THEN 

ENGCOST  =  .07833  *  KWHTOT(i) 

ELSEIF  KWHTOT(i)  >  180000  THEN 

ENGCOST  =  14099.4  +  ((KWHTOT(i)  -  180000)  *  .05247) 

ELSE 
END  IF 
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LASUMBILL: 

COSTl  =  MIN  +  DEMCOST  +  ENGCOST 
SURCOST  =  .05  *  COSTl 
BILL(i)  =  COSTl  +  SURCOST 
YRBILL  =  YRBILL  +  BILL(i) 
PMPENGY  =  PMPENGY  +  PMPKWH(i) 
FANENGY  =  FANENGY  +  FANKWH(i) 
CMPENGY  =  CMPENGY  +  CMPKWH(i) 
NEXT  i 


FOR  j  =  1  TO  12 

OPEN  ''B:LAEFF.DAT''  FOR  APPEND  AS  #4 
PRINT  #4,  "THE  MAX  MONTHLY  KW  IS 
"THE  BILLING  DEMAND 
"THE  MONTHLY  KWH 
"THE  MONTHLY  BILL 

•I  II 


PRINT  #4, 

PRINT  #4, 

PRINT  #4, 

PRINT  #4, 

PRINT  #4, 

IF  j  =  12  THEN 


=n 


=  tt 


KWPK(j) 
BILLDEM(j) 
KWHTOT( j) 
BILL(j) 


PRINT 

#4, 

"THE 

PRINT 

#4, 

"THE 

PRINT 

#4, 

"THE 

PRINT 

#4, 

"THE 

PRINT 

#4, 

"THE 

PRINT 

#4, 

II II 

PRINT 

#4, 

II II 

($)  =",  YRBILL 

WH)  =" ,  PMPENGY 
H)  FANENGY 


END  IF 
CLOSE  #4 
NEXT  j 
EXIT  SUB 

9 


NEWYORK: 

'It*************************************************************** 

•  **********************************************************ic***** 

'  CALCULATION  OF  MONTHLY  ENERGY  COSTS  IN  NEW  YORK  (CON.  ED.) 

'  Con  Edison  Rate  Schedule  No. 9 

'  General  Commercial  -  Large 

'  Eff.  Date  1  Apr.  1987 

'  Current  as  of  July  1989 

'  **  IN  NY  THE  MONTHLY  DEMAND  IS  STRICTLY  DET.  BY  THE  CURRENT 

'  MONTHS  PERFORMANCE.  THEREFORE  THE  DEMAND  COST  IS  AS  FOLLOWS; 

'***********************************************i,***i,i,*iti,i,i,i,i,i,ifi,if 

YRBILL  =  0 
PMPENGY  =  0 
FANENGY  =  0 
CMPENGY  =  0 
FOR  i  =  1  TO  12 

BILLDEM(i)  =  KWPK(i) 

DEMCOST  =  0 

IF  (i  >  5)  AND  (i  <  10)  GOTO  NYSUMMER 
MIN  =  74.18 
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IF  (BILLDEM(i)  >  5) 
DEMCOST  =  15.86  * 
ELSEIF  BILLDEM(i)  > 
DEMCOST  =  14274  + 
ELSE 
END  IF 

GOTO  NYENGYCOST 
NY SUMMER: 

MIN  =  96.98 
IF  (BILLDEM(i)  >  5) 
DEMCOST  =  20.36  * 
ELSEIF  BILLDEM(i)  > 
DEMCOST  =  18324  + 
ELSE 
END  IF 


AND  (BILLDEM{i)  <=  900)  THEN 
BILLDEM(i) 

900  THEN 

((BILLDEM(i)  -  900)  *  13.68) 


AND  (BILLDEM(i)  <=  900)  THEN 
BILLDEM(i) 

900  THEN 

((BILLDEM(i)  -  900)  *  18.18) 


>*********************************************************** 


'  CALCULATION  OF  THE  ENERGY  CHARGES 

f*************1t***1t*1fli1t**'k*************1t*******1c************ 


NYENGYCOST: 

IF  KWHTOT(i)  <=  15000  THEN 
ENGCOST  =  .0562  *  KWHTOT(i) 

ELSEIF  KWHTOT(i)  >  15000  THEN 

ENGCOST  =  843!  +  ((KWHTOT(i)  -  15000)  *  .0523) 
ELSE 
END  IF 


COSTl  =  MIN  +  DEMCOST  +  ENGCOST 
MAXBILL  =  .3207  *  KWHTOT(i) 

IF  COSTl  >  MAXBILL  THEN  COSTl  =  MAXBILL 
MINBILL  =  MIN 

IF  COSTl  <  MINBILL  THEN  COSTl  =  MINBILL 

/ 

FCADJ  =  .004414  *  KWHTOT(i) 

TAX  =  .0825  *  COSTl 

BILL(i)  =  COSTl  +  FCADJ  +  TAX 

YRBILL  =  YRBILL  +  BILL(i) 

PMPENGY  =  PMPENGY  +  PMPKWH(i) 

FANENGY  =  FANENGY  +  FANKWH(i) 

CMPENGY  =  CMPENGY  +  CMPKWH(i) 

NEXT  i 

f 

FOR  j  =  1  TO  12 


OPEN  *' 

'B: NYEFF.DAT”  FOR 

APPEND 

AS 

04 

PRINT 

04, 

"THE 

MAX  MONTHLY  KW 

IS 

S** 

9 

KWPK(j) 

PRINT 

04, 

"THE 

BILLING 

DEMAND 

s*’ 

9 

BILLDEM( j) 

PRINT 

04, 

"THE 

MONTHLY 

KWH 

^tt 

9 

KWHTOT(j) 

PRINT 

04, 

"THE 

MONTHLY 

BILL 

— 

9 

BILL(j) 

PRINT 

04, 

n  n 

PRINT 

04, 

fl  H 

IF  j  = 

■■  12 

THEN 

PRINT  04 

,  "THE  YEARLY 

TOTAL 

($) 

=",  YRBILL 
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IF  (BILLDEM(i)  >  5)  AND  (BILLDEM(i)  <=  900)  THEN 
DEMCOST  *  15.86  *  BILLDEM(i) 

ELSEIF  BILLDEM(i)  >  900  THEN 

DEMCOST  =  14274  +  ((BILLDEM(i)  -  900)  *  13.68) 
ELSE 
END  IF 

GOTO  NYENGYCOST 
NYSUMMER: 

MIN  =  96.98 

IF  (BILLDEM(i)  >  5)  AND  (BILLDEM(i)  <=  900)  THEN 
DEMCOST  =  20.36  *  BILLDEM(i) 

ELSEIF  BILLDEM(i)  >900  THEN 

DEMCOST  =  18324  +  ((BILLDEM(i)  -  900)  *  18.18) 
ELSE 
END  IF 

'  CALCULATION  OF  THE  ENERGY  CHARGES 
NYENGYCOST ; 

IF  KWHTOT(i)  <=  15000  THEN 
ENGCOST  =  .0562  *  KWHTOT(i) 

ELSEIF  KWHTOT(i)  >  15000  THEN 

ENGCOST  =  843!  +  ((KWHTOT(i)  -  15000)  *  .0523) 
ELSE 
END  IF 

/ 

COSTl  *  MIN  +  DEMCOST  +  ENGCOST 
MAXBILL  =  .3207  *  KWHTOT(i) 

IF  COSTl  >  MAXBILL  THEN  COSTl  =  MAXBILL 
MINBILL  =  MIN 

IF  COSTl  <  MINBILL  THEN  COSTl  =  MINBILL 

r 

FCADJ  =  .004414  *  KWHTOT(i) 

TAX  =  .0825  *  COSTl 

BILL(i)  =  COSTl  +  FCADJ  +  TAX 

YRBILL  =  YRBILL  +  BILL(i) 

PMPENGY  =  PMPENGY  +  PMPKWH(i) 

FANENGY  =  FANENGY  +  FANKWH{i) 

CMPENGY  =  CMPENGY  +  CMPKWH(i) 

NEXT  i 

FOR  j  =  1  TO  12 


OPEN  ' 

"B: NYEFF.DAT"  FOR 

APPEND 

AS 

#4 

PRINT 

#4, 

"THE 

MAX  MONTHLY  KW 

IS 

sit 

9 

KWPK(j) 

PRINT 

#4, 

"THE 

BILLING 

DEMAND 

9 

BILLDEM(j) 

PRINT 

#4, 

"THE 

MONTHLY 

KWH 

9 

KWHTOT(j) 

PRINT 

#4, 

"THE 

MONTHLY 

BILL 

9 

BILL(-j) 

PRINT 

#4, 

nn 

PRINT 

#4, 

nit 

IF  j  = 

=  12 

THEN 

PRINT  #4,  "THE  YEARLY  TOTAL  ($)  YRBILL 
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SELECT  CASE  PMPHP 
CASE  0  TO  5 
PMPCST  =  2690 
VFDCST  =  5183 
CASE  IS  <=  7.5 
PMPCST  =  2740 
VFDCST  =  5183 
CASE  IS  <=  10 
PMPCST  =  3275 
VFDCST  =  5183 
CASE  IS  <=  15 
PMPCST  =  3765 
VFDCST  =  6933 
CASE  IS  <=  20 
PMPCST  =  4365 
VFDCST  =  8618 
CASE  IS  <=  25 
PMPCST  =  4440 
VFDCST  =  10738 
CASE  IS  <=  30 
PMPCST  =  5365 
VFDCST  =  11938 
CASE  IS  <=  50 
PMPCST  =  6795 
VFDCST  =  13948 
CASE  IS  >  50 
PMPCST  =  8720 
VFDCST  =  18050 

END  SELECT 

/ **************************************************************** 

'  DETERMINE  THE  COST  OF  THE  ELECT  INSTALLATION  FOR  THE  TOWER  FAN 

'**************************************************************** 

SELECT  CASE  FANHP 
CASE  0  TO  5 
FANCST  =  1315 
CASE  IS  <=  7.5 
FANCST  =  1315 
CASE  IS  <=  10 
FANCST  =  1450 
CASE  IS  <=  15 
FANCST  =  1815 
CASE  IS  <=  20 
FANCST  =  1965 
CASE  IS  <=  25 
FANCST  =  1965 
CASE  IS  <=  30 
FANCST  =  2500 
CASE  IS  <=  50 
FANCST  =  2750 
CASE  IS  >  50 
FANCST  =  4085 

END  SELECT 
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•  *********************************************************** 

'  DETERMINE  THE  CAPITAL  COST  OF  THE  PIPING , VALVING , FITTINGS 
•  *********************************************************** 
SELECT  CASE  diameter 
CASE  5.05 
PIPECST  =  19850 
CASE  6.07 
PIPECST  =  23495 
CASE  8.07 
PIPECST  =  33548 
CASE  10.02 
PIPECST  =  52217 
END  SELECT 

'  ********************************************************** 

'  DETERMINE  THE  TOTAL  CAPITAL  COST 

I  ********************************************************** 

TOTCST  =  TWRCST  +  PMPCST  +  FANCST  +  PIPECST  +  VFDCST 

OVRHD  =  .15  *  TOTCST 

SUBl  =  TOTCST  +  OVRHD 

PROF  =  .1  *  SUBl 

CAPCST  =  SUBl  +  PROF 

•  ********************************************************** 

•  INCLUDE  THE  1989  MEANS  CITY  COST  INDEXES 

•  ********************************************************** 

SELECT  CASE  CITY 
CASE  1 

CTYINDX  =  .917 
CASE  2 

CTYINDX  =  1.117 
CASE  3 

CTYINDX  =  1.235 
END  SELECT 

CAPCST  =  CAPCST  *  CTYINDX 
SELECT  CASE  CITY 
CASE  1 

OPEN  "BtATLEFF.DAT"  FOR  APPEND  AS  #4 
CASE  2 

OPEN  "B: LAEFF.DAT"  FOR  APPEND  AS  #4 
CASE  3 

OPEN  "B: NYEFF.DAT"  FOR  APPEND  AS  #4 
END  SELECT 


PRINT 

#4, 

"THE 

TOWER  COST  IS 

ff 

9 

TWRCST 

PRINT 

"THE 

PUMP  COST  IS 

It 

9 

PMPCST 

PRINT 

#4, 

"THE 

FAN  ELECT  COST", 

FANCST 

PRINT 

#4, 

"THE 

PIPE  SYS  COST 

If 

r 

PIPECST 

PRINT 

#4, 

"THE 

SUBTOTAL  IS 

ft 

9 

TOTCST 

PRINT 

#4, 

"THE 

OVERHEAD 

H 

9 

OVRHD 

PRINT 

#4, 

"THE 

PROFIT  IS 

tf 

9 

PROF 

PRINT 

#4, 

"THE 

CITY  INDEX  IS 

If 

9 

CTYINDX 

PRINT 

/4, 

"THE 

TOTAL  COST  IS 

ff 

9 

CAPCST 

CLOSE 

#4 

END  SUB 
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APPENDIX  B 


DATA  GENERATED  FROM  WEBER  MODEL  COMPARISON 
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OUTPUT  OF  DELTA  TEMPERATURE  FOR  EFFECTIVENESS  AND  JOYCE  METHOD 


DELTA  EFFECTIVENESS 

DELTA  WEBER 

DELTA  WEBER 

DELTA  WEBER 

(DEG  F) 

(DEG  F) 

PLUS  15% 

MINUS  15% 

1.6 

1.4 

1.61 

1.19 

6.4 

5.9 

6.79 

5.02 

6.5 

5.9 

6.79 

5.02 

6.5 

5.9 

6.79 

5.02 

7.8 

7.3 

8.39 

6.2 

9.8 

8.9 

10.235 

7.57 

11 

10 

11.5 

8.5 

11 

10 

11.5 

8.5 

11.1 

10.4 

11.96 

8.84 

12.2 

11.2 

12.88 

9.52 

14.5 

13 

14.95 

11.05 

14.6 

13 

14.95 

11.05 

14.6 

13 

14.95 

11.05 

19.9 

17.8 

20.47 

15.13 

20 

18.1 

20.82 

15.39 

20.1 

17.8 

20.47 

15.13 

20.2 

17.8 

20.47 

15.13 

26.9 

23.4 

26.91 

19.89 

26.9 

23.4 

26.91 

19.89 

26.9 

23.4 

26.91 

19.89 

27 

23.6 

27.14 

20.06 

28.6 

25 

28.75 

21.25 

28.8 

25 

28.75 

21.25 

29 

25 

28.75 

21.25 
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APPENDIX  C 


DATA  GENERATED  FROM  JOYCE  MODEL  COMPARISON 
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AT750.DAT  PDNP  SIZE  750  6PM  JOYCE  METHOD 


MONTH=  12 

BIN  *  57 

START  TIME 

=  11:12:25 

END  TIME  =  11:53 

:  09  RUNTIME  = 

40:44 

Ctrl=  7 

Ctr2  =  92 

Ctr3  =  3 

1=  12 

J=  10 

TDB(DEG  F) 

TWB(DEG 

F)  HOURS 

TOWER 

57 

52 

24.82143 

3 

GPM 

%GPM 

CFM 

TONS 

750 

100 

120000 

100 

CWS(DEG  F) 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER 

TOWERDUTY 

80 

84.03179 

4.031794 

19.60413 

.1174005 

PIPE  D(FT) 

Pipevelocity (FT/SEC)  REd 

pipe  frict 

.6666667 

4.787051 

369371.3 

1.592428E-02 

pipeL  total  PHpipe 

PHcond 

PHtower 

PHtotal(FT 

Water) 

666.6667 

5.666439 

5.105591 

11 

21.77203 

PASSES 

Condenserve loc i ty 

2 

2.900009 

UACOND 

LMTDCOND 

NfISENT 

NFTON 

COND.TEMP 

893100.2 

1.690421 

.5164742 

.6406521 

84.4407 

KWTONGUESS 

KWTONMODEL 

KWHBIN 

KWHTOT 

.9074575 

.9074562 

2526.094 

6573.959 

HOURS 

HOURSSUM 

24.82143 

57.5 

CHILLERKW 

PUMPKW 

FANKW 

TOTALKW 

90.74575 

6.151319 

4.87361 

101.7707 

THE  YEARLY  TOTAL  ($)  =  334021.2 

THE  PUMP  ENERGY  (KWH)  =  17345.86 

THE  FAN  ENERGY  (KWH)  =  46195.46 

THE  CHILLER  ENERGY  (KWH)  =  337459.7 

THE  TOTALSYSTEM  ENERGY  (KWH)  =  401001 
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AT100.DAT  PUMP  SIZE  1000  6PM  JOYCE  METHOD 


MONTH=  12  BIN  =  57 

START  TIME  =  11:55:00  END  TIME  =  12:37:03  RUNTIME  =  42:03 

ctrl=  7  ctr2  =  92  ctr3  =  3 


1=  12 

TDB(DEG  F) 
57 
GPM 
750 

CWS(DEG  F) 
80 

PIPE  D(FT) 
.6666667 


J=  10 

TWB(DEG  F) 
52 

%GPM 

75 

CWR(DEG  F) 
84.03179 


HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 


Pipevelocity (FT/SEC)  REd 
4.787051  369371.3 


plpeL  total  PHpipe 
PHtotal(FT  Water) 


PHcond 


666.6667 

PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074575 

HOURS 

24.82143 

CHILLERKW 

90.74575 


5.666439  5.105591 

Con den s er ve 1 oc i ty 
2.900009 


LMTDCOND 

1.690421 

KWTONMODEL 

.9074562 

HOURSSUM 

57.5 

PUMPKW 

4.342596 


NfISENT 

.5164742 

KWHBIN 

2481.199 


FANKW 

4.87361 


TOWER 

3 

TONS 

100 

DELTATTOWER  TOWERDUTY 

19.60413  .1174005 

pipe  frict 
1.592428E-02 
PHtower 


11  21.77203 


NFTON  COND . TEMP 

.6406521  84.4407 

KWHTOT 
6469.957 


TOTALKW 

99.96196 


THE  YEARLY  TOTAL  ($)  = 

THE  PUMP  ENERGY  (KWH)  » 

THE  FAN  ENERGY  (KWH)  * 

THE  CHILLER  ENERGY  (KWH)  * 
THE  TOTALSYSTEM  ENERGY  (KWH) 


333374 

9823.913 

46195.46 

337459.7 


393479.1 


43 


AT125.DAT  PUMP  8I2E  1250  OPM  JOYCE  METHOD 


MONTH=  12 
START  TIME 
Ctrl=  7 
1=  12 

TDB(DEG  F) 
57 
GPM 
750 

CWS(DEG  F) 
80 

PIPE  D(FT) 
.6666667 


BIN  *  57 

11:36:25  END  TIME  =  12:19:33  RUNTIME  =  43:08 
ctr2  *  92  Ctr3  =  3 


J=  10 

TWB(DEG  F) 
52 

%GPM 

60 

CWR(DEG  F) 
84.03179 


HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 


Pipevelocity (FT/SEC)  REd 
4.787051  369371.3 


TOWER 

3 

TONS 

100 

DELTATTOWER  TOWERDUTY 

19.60413  .1174005 

pipe  frict 
1.592428E-02 


PHpipe  PHcond 

Water) 

5.666439  5.105591 

Condenservelocity 


pipeL  total 
PH total (FT 
666.6667 
PASSES 
2 

UACOND 

893100.2 

KWTONGUESS 

.9074575 

HOURS 

24.82143 

CHILLERKW 

90.74575 


2.900009 

LMTDCOND 

1.690421 

KWTONMODEL 

.9074562 

HOURSSUM 

57.5 

PUMPKW 

4.541853 


NfISENT 

.5164742 

KWHBIN 

2486.145 


FANKW 

4.87361 


PHtower 


11 


NFTON 

.6406521 

KWHTOT 

6481.415 


TOTALKW 

100.1612 


21.77203 

COND.TEMP 

84.4407 


THE  YEARLY  TOTAL  ($)  = 

THE  PUMP  ENERGY  (KWH)  = 

THE  F7VN  ENERGY  (KWH)  = 

THE  CHILLER  ENERGY  (KWH)  = 
THE  TOTALSYSTEM  ENERGY  (KWH) 


333380.8 

10126.37 

46195.46 

337459.7 


393781.5 
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AT150.DAT  PDMP 

8I2E  1500  GPM 

JOYCE  METHOD 

MONTH*  12 

BIN  =  57 

START  TIME 

*  12:25:32  END  TIME  *  13:09 

:23  RUNTIME  * 

43:51 

Ctrl*  7 

ctr2  *  92 

ctr3  *  3 

I*  12 

J*  10 

TDB(DEG  F) 

TWB(DEG  F) 

HOURS 

TOWER 

57 

52 

24.82143 

3 

GPM 

%GPM 

CFM 

TONS 

750 

50 

120000 

100 

CWS(DEG  F) 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER 

TOWERDUTY 

80 

84.03179 

4.031794 

19.60413 

.1174005 

PIPE  D(FT) 

Pipevelocity (FT/SEC)  REd 

pipe  frict 

.6666667 

4.787051 

369371.3 

1.592428E-02 

pipeL  total  PHpipe 

PHcond 

PHtower 

PHtotal(FT 

Water) 

21.77203 

666.6667 

5.666439 

5.105591 

11 

PASSES 

Condenservelocity 

2 

2.900009 

UACOND 

LMTDCOND 

NfISENT 

NFTON 

COND.TEMP 

893100.2 

1.690421 

.5164742 

.6406521 

84.4407 

KWTONGUESS 

KWTONMODEL 

KWHBIN 

KWHTOT 

* 

.9074575 

.9074562 

2501.478 

6516.936 

HOURS 

HOURSSUM 

24.82143 

57.5 

CHILLERKW 

PUMPKW 

FANKW 

TOTALKW 

90.74575 

5.159607 

4.87361 

100.779 

THE  YEARLY 

TOTAL  ($)  = 

333485.1 

THE  PUMP  ENERGY  (KWH)  = 

11436.89 

THE  FAN  ENERGY  (KWH)  = 

46195.46 

THE  CHILLER  ENERGY  (KWH) 

*  337459.7 

THE  TOTALSYSTEM  ENERGY  (KWH)  * 

395092 
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ATEF750.DAT  PUMP  SIZE  750  6PM  EFFECTZVEME8S  METHOD 


MONTH*  12  BIN  *  57 

TIME  START  =  16:44:50  TIME  END  =  16:51:52  RUNTIME  =7:02 
ctrl=  7  ctr2  =  92  Ctrl  =  3 


1=  12 

TDB(DEG  F) 

57 

UPM 

750 

CWS(DEG  F) 

80 

9.490044E-02 
PIPE  D(FT) 

. 6666667 


J=  10 

TWB(DEG  F) 

HOURS 

TOWER 

52 

24.82143 

3 

%GPM 

CFM 

TONS 

100 

120000 

100 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER  TOWERDUTY 

84.03179 

4.031794  21.74524 

Pipevelocity (FT/SEC)  REd 

pipe  frict 

4.787051 

369371.3 

1.592428E-02 

pipeL  total  PHpipe 
PHtotal(FT  Water) 


PHcond 


666.6667 

PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074576 

HOURS 

24.82143 

CHILLERKW 

90.74576 


5.666439  5.105591 

Condenservelocity 
2.900009 


LMTDCOND 

1.690421 

KWTONMODEL 

.9074562 

HOURSSUM 

57.5 

PUMPKW 

6.151319 


NfISENT 

.5164742 

KWHBIN 

2502.91 


FANKW 

3.939574 


PHtOwer 


11  21.77203 


NFTON  COND . TEMP 

.6406521  84.4407 

KWHTOT 
6502.745 


TOTALKW 

100.8366 


THE  YEARLY  TOTAL  ($)  = 

THE  PUMP  ENERGY  (KWH)  = 

THE  FAN  ENERGY  (KWH)  = 

THE  CHILLER  ENERGY  (KWH)  = 
THE  TOTALSYSTEM  ENERGY  (KWH) 


333296 

17345.86 

40526.96 

337296.7 


395169.5 
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ATEFIOO 

.DAT  PUMP  SIZE 

1000  GPM 

EFFECTIVENESS  METHOD 

MONTH=  12 

BIN  =  57 

TIME  START  = 

17:22:58  TIME  END  =  17:29:55  RUNTIME  = 

6:57 

ctrl=  7 

ctr2  =  92 

ctr3  *  3 

1=  12 

J=  10 

TDB{DEG  F) 

TWB(DEG  F) 

HOURS 

TOWER 

57 

52 

24.82143 

3 

GPM 

%GPM 

CFM 

TONS 

750 

75 

120000 

100 

CWS(DEG  F) 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER 

TOWERDUTY 

80 

84.03179 

4.031794  21.74524 

9.490044E-02 

PIPE  D(FT) 

Pipevelocity (FT/SEC)  REd 

pipe  frict 

.6666667 

4.787051 

369371.3 

1.592428E-02 

pipeL  total 

PHpipe 

PHcond 

PHtower 

PHtotal(FT  Water) 

666.6667 

5.666439 

5.105591 

11 

21.77203 

PASSES 

Condenservelocity 

2 

2.900009 

UACOND 

LMTDCOND 

NfISENT 

NFTON 

COND . TEMP 

893100.2 

1.690421 

.5164742 

.6406521 

84.4407 

KWTONGUESS 

KWTONMODEL 

KWHBIN 

KWHTOT 

.9074576 

.9074562 

2458.015 

6398.744 

HOURS 

HOURSSUM 

24.82143 

57.5 

CHILLERKW 

PUMPKW 

FANKW 

TOTALKW 

90.74576 

4.342596 

3.939574 

99.02793 

THE  YEARLY  TOTAL  ($)  = 

332648.8 

THE  PUMP  ENERGY  (KWH)  = 

9823.913 

THE  FAN  ENERGY  (KWH)  = 

40526.96 

THE  CHILLER  ENERGY  (KWH)  = 

337296.7 

THE  TOTALS YSTEM  ENERGY  (KWH) 

387647.6 
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ATEF125.DAT  PUMP  SIZE  1250  6PM  EFFECTIVEMEbS  METHOD 


MONTH*  12 

BIN  =  57 

TIME  START  = 

17:32:44  TIME  END  *  17:40: 

57  RUNTIME  =8:13 

Ctrl*  7 

Ctr2  =  92 

ctr3  =  3 

1=  12 

J=  10 

TDB(DEG  F) 

TWB(DEG  F) 

HOURS 

TOWER 

57 

52 

24.82143 

3 

GPM 

%GPM 

CFM 

TONS 

750 

60 

120000 

100 

CWS(DEG  F) 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER  TOWERDUTY 

80 

84.03179 

4.031794  21.74524 

9.490044E-02 

PIPE  D(FT) 

Pipevelocity (FT/SEC)  REd 

pipe  frict 

.6666667 

4.787051 

369371.3 

1.592428E-02 

pipeL  total 

PHpipe 

PHcond 

PHtower 

PHtotal(FT  Water) 

666.6667 

5.666439 

5.105591 

11  21.77203 

PASSES 

Condenservelocity 

2 

2.900009 

UACOND 

LMTDCOND 

NfISENT 

NFTON  COND . TEMP 

893100.2 

1.690421 

.5164742 

.6406521  84.4407 

KWTONGUE3S 

KWTONMODEL 

KWHBIN 

KWHTOT 

.9074576 

.9074562 

2462.961 

6410.201 

HOURS 

HOURSSUM 

24.82143 

57.5 

CHILLERKW 

PUMPKW 

FANKW 

TOTALKW 

90.74576 

4.541853 

3.939574 

99.22719 

THE  YEARLY  TOTAL  ($)  =  332655.6 

THE  PUMP  ENERGY  (KWH)  =  10126.37 

THE  FAN  ENERGY  (KWH)  =  40526.96 

THE  CHILLER  ENERGY  (KWH)  =  337296.7 

THE  TOTALSYSTEM  ENERGY  (KWH)  =  387950 
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ATEF150.DAT  PDMP  SIZE  1500  GPM  EFFECTIVENESS  METHOD 


MONTH=  12  BIN  =  57 

TIME  START  =  17:42:20  TIME  END  =  17:51:43  RUNTIME  *  9:23 
ctrl=  7  Ctr2  =  92  Ctr3  =  3 


1=  12 

TDB(DEG  F) 

57 

GPM 

750 

CWS(DEG  F) 

80 

9.490044E-02 
PIPE  D(FT) 
.6666667 


J=  10 

TWB(DEG  F) 

HOURS 

TOWER 

52 

24.82143 

3 

%GPM 

CFM 

TONS 

50 

120000 

100 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER  TOWERDUTY 

84.03179 

4.031794  21.74524 

Pipevelocity (FT/SEC)  REd 

pipe  frict 

4.787051 

369371.3 

1.592428E-02 

pipeL  total  PHpipe 
PHtotal(FT  Water) 
666.6667  5.666439 


PHcond 


PHtower 


5.105591  11 


21.77203 


Condenservelocity 


PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074576 

HOURS 

24.82143 

CHILLERKW 

90.74576 


2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
5.159607 


NfISENT 

.5164742 

KWHBIN 

2478.294 


FANKW 

3.939574 


NFTON  COND . TEMP 

.6406521  84.4407 

KWHTOT 
6445.722 


TOTALKW 

99.84494 


THE  YEARLY  TOTAL  ($)  = 

THE  PUMP  ENERGY  (KWH)  = 

THE  FAN  ENERGY  (KWH)  = 

THE  CHILLER  ENERGY  (KWH)  = 
THE  TOTALSYSTEM  ENERGY  (KWH) 


332759.9 

11436.89 

40526.96 

337296.7 


389260.5 
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ATT0W2.DAT  TOWER  #  2  JOYCE  METHOD 


MONTH=  12 

BIN  =  57 

START  TIME 

=  10:02:55 

END  TIME  =10 

:43:28  RUNTIME  =  40:33 

ctrl=  7 

ctr2  =  90 

Ctr3  =  3 

1=  12 

J=  10 

TDB(DEG  F) 

TWB(DEG  F) 

HOURS 

TOWER 

57 

52 

24.82143 

2 

GPM 

%GPM 

CFM 

TONS 

750 

75 

120000 

100 

CWS(DEG  F) 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER 

TOWERDUTY 

80 

84.03179 

4.031794 

19.00449 

.1246106 

PIPE  D(FT) 

Pipevelocity(FT/SEC)  REd 

pipe  frict 

.6666667 

4.787051 

369371.3 

1.592428E-02 

pipeL  total  PHpipe 

PHcond 

PHtower 

PHtOtal(FT 

Water) 

666.6667 

5.666439 

5.105591 

11 

21.77203 

PASSES 

Condenservelocity 

2 

2.900009 

UACOND 

LMTDCOND 

NfISENT 

NFTON 

COND . TEMP 

893100.2 

1.690421 

.5164742 

.6406521 

84.4407 

KWTONGUESS 

KWTONMODEL 

KWHBIN 

KWHTOT 

.9074574 

.9074562 

2541.575 

6686.729 

HOURS 

HOURSSUM 

24.82143 

57.5 

CHILLERKW 

PUMPKW 

FANKW 

TOTALKW 

90.74574 

4.342596 

7.306037 

102.3944 

THE  YEARLY  TOTAL  ($)  = 

THE  PUMP  ENERGY  (KWH)  = 

THE  FAN  ENERGY  (KWH)  = 

THE  CHILLER  ENERGY  (KWH)  = 
THE  TOTALSYSTEM  ENERGY  (KWH) 


336101.8 
9823.913 
67227.36 

337547.8 


414599.1 


ATT0W5.DAT  TOWER  #  5  JOYCE  METHOD 


MONTH=  12 
START  TIME  = 
ctrl=  7 
1=  12 

TDB(DEG  F) 
57 
GPM 
750 

CWS(DEG  F) 
80 

PIPE  D(FT) 
.6666667 
pipeL  total 
PHtotal(FT 
666.6667 
PASSES 
2 

UACOND 

893100.2 

KWTONGUESS 

.9074575 

HOURS 

24.82143 

CHILLERKW 

90.74575 


BIN  =  57 
10:44:54 
Ctr2  =  92 

J=  10 

TWB(DEG  F) 
52 

%GPM 

75 

CWR(DEG  F) 
84.03179 


2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURSSUM 
57.5 
PUMPKW 
4.342596 


END  TIME  =  11 
Ctr3  =  3 

HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 


NfISENT 

.5164742 

KWHBIN 

2481.199 


FANKW 

4.87361 


Pipevelocity (FT/SEC)  REd 
4.787051  369371.3 

PHpipe  PHcond 

Water) 

5.666439  5.105591 

Condenservelocity 


26:56  RUNTIME  =  42:02 


TOWER 

5 

TONS 

100 

DELTATTOWER  TOWERDUTY 
19.60413  .1174005 

pipe  frict 
1.592428E-02 
PHtower 


11 


NFTON 

.6406521 

KWHTOT 

6469.957 


TOTALKW 

99.96196 


21.77203 

COND.TEMP 

84.4407 


THE  YEARLY  TOTAL  ($)  = 

THE  PUMP  ENERGY  (KWH)  = 

THE  FAN  ENERGY  (KWH)  = 

THE  CHILLER  ENERGY  (KWH)  = 
THE  TOTALS Y STEM  ENERGY  (KWH) 


333374 

9823.913 

46195.46 

337459.7 


393479.1 
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ATTOir9.DAT  TOWER  #  9 


JOYCE  METHOD 


MONTH=  12 
START  TIME 
Ctrl=  7 
1=  12 

TDB(DEG  F) 
57 
GPM 
750 

CWS(DEG  F) 
80 

PIPE  D(FT) 
.6666667 


BIN  =  57 

11:29:52  END  TIME  =  12:13:13  RUNTIME  =  43:21 
ctr2  =  91  ctr3  =  3 


J*  10 

TWB(DEG  F) 
52 

%GPM 

75 

CWR(DEG  F) 
84.03179 


HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 


Pipevelocity (FT/SEC)  REd 
4.787051  369371.3 


TOWER 

9 

TONS 

100 

DELTATTOWER  TOWERDUTY 
19.93193  .1136424 

pipe  frict 
1.592428E-02 


pipeL  total  PHpipe  PHcond 

PHtotal(FT  Water) 

666.6667  5.666439  5.105591 

PASSES  Condenservelocity 


PHtower 


2 

UACOND 

893100.2 

KWTONGUESS 

.9074575 

HOURS 

24.82143 

CHILLERKW 

90.74575 


2.900009 

LMTDCOND 

1.690421 

KWTONMODEL 

.9074562 

HOURSSUM 

57.5 

PUMPKW 

4.342596 


NfISENT 

.5164742 

KWHBIN 

2405.123 


FANKW 

1.808699 


11 


NFTON 

.6406521 

KWHTOT 

6187.521 


TOTALKW 

96.89704 


21.77203 

COND.TEMP 

84.4407 


THE  YEARLY  TOTAL  ($)  = 

THE  PUMP  ENERGY  (KWH)  = 

THE  FAN  ENERGY  (KWH)  = 

THE  CHILLER  ENERGY  (KWH)  = 
THE  TOTALSYSTEM  ENERGY  (KWH) 


329604.7 

9823.913 

17438.78 

337425.3 


364687.9 
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ATT0W12.DAT  TOWER  #  12  JOYCE  METHOD 


MONTH=  12 
START  TIME 
Ctrl=  7 
1=  12 

TDB(DEG  F) 
57 
GPM 
750 

CWS(DEG  F) 
80 

PIPE  D(FT) 
.6666667 


BIN  =  57 

12:21:59  END  TIME  =  13:02:36  RUNTIME  =  40:37 
Ctr2  =  91  Ctr3  *  3 


J=  10 

TWB(DEG  F) 
52 

%GPM 

75 

CWR(DEG  F) 
84.03179 


HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 


Pipevelocity (FT/SEC)  REd 
4.787051  369371.3 


TOWER 

12 

TONS 

100 

DELTATTOWER  TOWERDUTY 
19.7356  .1158782 

pipe  frict 
1.592428E-02 


pipeL  total  PHpipe  PHcond 

PHtotal(FT  Water) 

666.6667  5.666439  5.105591 

PASSES  Condenservelocity 


PHtower 


2 

UACOND 

893100.2 

KWTONGUESS 

.9074575 

HOURS 

24.82143 

CHILLERKW 

90.74575 


2.900009 

LMTDCOND 

1.690421 

KWTONMODEL 

.9074562 

HOURSSUM 

57.5 

PUMPKW 

4.342596 


NfISENT 

.5164742 

KWHBIN 

2397.792 


FANKW 

1.513313 


11 


NFTON 

.6406521 

KWHTOT 

6x59.366 


TOTALKW 

96.60166 


21.77203 

COND.TEMP 

84.4407 


THE  YEARLY  TOTAL  ($)  = 

THE  PUMP  ENERGY  (KWH)  = 

THE  FAN  ENERGY  (KWH)  » 

THE  CHILLER  ENERGY  (KWH)  = 
THE  TOTALS YSTEM  ENERGY  (KWH) 


329208.7 

9823.913 

14440.33 

337443.6 


361707.9 
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ATT0W16.DAT 


TOWER  #  16  JOYCE  METHOD 


MONTH=  12  BIN  =  57 

START  TIME  =  13:04:09  END  TIME  =  13:46:03  RUNTIME  =  41:54 

Ctrl-  7  ctr2  =  90  ctr3  *  3 


1=  12 

TDB(DEG  F) 
57 
GPM 
750 

CWS(DEG  F) 
80 

PIPE  D(FT) 
.6666667 
pipeL  total 


HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 


J=  10 

TWB(DEG  F) 

52 

%GPM 
75 

CWR(DEG  F) 

84.03179 

Pipevelocity (FT/SEC)  REd 
4.787051  369371.3 

PHplpe  PHcond 

PHtotal(FT  Water) 

666.6667  5.666439  5.105591 

Condenservelocity 
2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURSSUM 
57.5 
PUMPKW 
4.342596 


PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074575 

HOURS 

24.82143 

CHILLERKW 

90.74575 


NfISENT 

.5164742 

KWHBIN 

2385.629 


FANKW 

1.023324 


TOWER 

16 

TONS 

100 

DELTATTOWER  TOWERDUTY 
20.19933  .1106671 

pipe  frict 
1.592428E-02 
PHtower 


11  21.77203 


NFTON  COND . TEMP 

.6406521  84.4407 

KWHTOT 
6114.878 


TOTALKW 

96.11167 


THE  YEARLY  TOTAL  ($)  = 

THE  PUMP  ENERGY  (KWH)  = 

THE  FAN  ENERGY  (KWH)  = 

THE  CHILLER  ENERGY  (KWH)  = 
THE  TOTALSYSTEM  ENERGY  (KWH) 


328628.1 
9823.913 
10003.45 

337395.2 


357222.5 
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ATTOW19.DAT 

TOWER  #  19 

JOYCE  METHOD 

MONTH*  12 

BIN  *  57 

START  TIME 

=  14:35:43 

END  TIME  *  15: 

19:20  RUNTIME 

=  43:37 

Ctrl*  7 

Ctr2  *  90 

ctr3  *  3 

1=  12 

J=  10 

TDB(DEG  F) 

TWB(DEG  F) 

HOURS 

TOWER 

57 

52 

24.82143 

19 

GPM 

%GPM 

CFM 

TONS 

750 

75 

120000 

100 

CWS(DEG  F) 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER 

TOWERDUTY 

80 

84.03179 

4.031794 

20.77497 

.104522 

PIPE  D(FT) 

Pipevelocity (FT/SEC)  REd 

pipe  frict 

.6666667 

4.787051 

369371.3 

1.592428E-02 

pipeL  total  PHpipe 

PHcond 

PHtower 

PHtotal(FT 

Water) 

666.6667 

5.666439 

5.105591 

11 

21.77203 

PASSES 

Condenservelocity 

2 

2.900009 

UACOND 

LMTDCOND 

NfISENT 

NFTON 

COND.TEMP 

893100.2 

1.690421 

.5164742 

.6406521 

84.4407 

KWTONGUESS 

KWTONMODEL 

KWHBIN 

KWHTOT 

.9074577 

.9074562 

2385.587 

6115.9 

HOURS 

HOURSSUM 

24.82143 

57.5 

CHILLERKW 

PUMPKW 

FANKW 

TOTALKW 

90.74577 

4.342596 

1.021625 

96.10999 

THE  YEARLY 

TOTAL  ($)  = 

328658.3 

THE  PUMP  ENERGY  (KWH)  = 

9823.913 

THE  FAN  ENERGY  (KWH)  * 

10252.19 

THE  CHILLER  ENERGY  (KWH) 

=  337320.8 

THE  TOTALS Y STEM  ENERGY  (KWH)  * 

357396.9 
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ATEF2.DAT  TOWER  #  2  EFFECTIVENESS  METHOD 


MONTH*  12 

BIN  *  57 

TIME  START 

*  15:37:19 

TIME  END  *  15:44 

:  00  RUNTIME  = 

=  6:41 

Ctrl*  7 

Ctr2  =  91 

ctr3  *  3 

I*  12 

J*  10 

TDB(DEG  F) 

TWB(DEG 

F)  HOURS 

TOWER 

57 

52 

24.82143 

2 

GPM 

%GPM 

CFM 

TONS 

750 

75 

120000 

100 

CWS(DEG  F) 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER 

TOWERDUTY 

80 

84.03179 

4.031794 

20.96244 

.1025935 

PIPE  D(FT) 

Pipevelocity (FT/SEC)  REd 

pipe 

frict 

.6666667 

4.787051 

369371.3 

1.592428E- 

•02 

plpeL  total  PHpipe 

PHcond 

PHtower 

PHtotal(FT 

Water) 

666.6667 

5.666439 

5.105591 

11 

21.77203 

PASSES 

Condenservelocity 

2 

2.900009 

UACOND 

LMTDCOND 

NfISENT 

NFTON 

COND . TEMP 

893100.2 

1.690421 

.5164742 

.6406521 

84.4407 

KWTONGUESS 

KWTONMODEL 

KWHBIN 

KWHTOT 

.9074575 

.9074562 

2509.533 

6588.154 

HOURS 

HOURS SUM 

24.82143 

57.5 

CHILLERKW 

PUMPKW 

FANKW 

TOTALKW 

90.74575 

4.342596 

6.015156 

101.1035 

THE  YEARLY  TOTAL  ($)  =  335101.8 

THE  PUMP  ENERGY  (KWH)  =  9823.913 

THE  FAN  ENERGY  (KWH)  =  59464.87 

THE  CHILLER  ENERGY  (KWH)  =  337298.5 

THE  TOTALSYSTEM  ENERGY  (KWH)  =  406587.3 
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ATEF5.DAT  TOWER  #  5  EFFECTIVENESS  METHOD 


MONTH=  12  BIN  =  57 

TIME  START  =  15:46:17  TIME  END  *  15:54:27  RUNTIME  =8:10 
ctrl=  7  ctr2  =  92  Ctrl  =  3 


1=  12 

TDB(DEG  F) 

57 

GPM 

750 

CWS(DEG  F) 

80 

9.490044E-02 
PIPE  D(FT) 
.6666667 


J=  10 

TWB(DEG  F) 

HOURS 

TOWER 

52 

24.82143 

5 

%GPM 

CFM 

TONS 

75 

120000 

100 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER  TOWERDUTY 

84.03179 

4.031794  21.74524 

Pipevelocity(FT/SEC)  REd 

pipe  frict 

4.787051 

369371.3 

1.592428E-02 

pipeL  total  PHpipe  PHcond 

PHtotal(FT  Water) 

666.6667  5.666439  5.105591 

PASSES  Condenservelocity 


PHtower 


2 

UACOND 

893100.2 

KWTONGUESS 

.9074576 

HOURS 

24.82143 

CHILLERKW 

90.74576 


2.900009 

LMTDCOND 

1.690421 

KWTONMODEL 

.9074562 

HOURSSUM 

57.5 

PUMPKW 

4.342596 


NfISENT 

.5164742 

KWHBIN 

2458.015 


FANKW 

3.939574 


11 


NFTON 

.6406521 

KWHTOT 

6398.744 


TOTALKW 

99.02793 


21.77203 

COND.TEMP 

84.4407 


THE  YEARLY  TOTAL  ($)  = 

THE  PUMP  ENERGY  (KWH)  = 

THE  FAN  ENERGY  (KWH)  = 

THE  CHILLER  ENERGY  (KWH)  = 
THE  TOTALSYSTEM  ENERGY  (KWH) 


332648.8 

9823.913 

40526.96 

337296.7 


387647.6 


57 


ATEF9.DAT  TOWER  #  9  EFFECTIVENESS  METHOD 


MONTH=  12  BIN  =  57 

TIME  START  =  15:57:15  TIME  END  =  16:06:35  RUNTIME  =9:20 
ctrl=  7  ctr2  =  92  ctr3  =  3 

1=  12  J=  10 

TDB(DEG  F)  TWB(DEG  F)  HOURS  TOWER 

57  52  24.82143  9 

GPM  %GPM  CFM  TONS 

750  75  120000  100 

CWS(DEG  F)  CWR(DEG  F)  DELTATCOND  DELTATTOWER  TOWERDUTY 

80  84.03179  4.031794  22.17689 

9.089065E-02 

PIPE  D(FT)  Pipevelocity (FT/SEC)  REd  pipe  frict 

.6666667  4.787051  369371.3  1.592428E-02 

pipeL  total  PHpipe  PHcond  PHtower 

PHtotal(FT  Water) 

666.6667  5.666439  5.105591  11  21.77203 

PASSES  Condenservelocity 

2  2.900009 

UACOND  LMTDCOND  NfISENT  NFTON  COND.TEMP 

893100.2  1,690421  .5164742  .6406521  84.4407 

KWTONGUESS  KWTONMODEL  KWHBIN  KWHTOT 

.9074577  .9074562  2396.136  6159.948 

HOURS  HOURSSUM 

24.82143  57.5 

CHILLERKW  PUMPKW  FANKW  TOTALKW 

90.74577  4.342596  1.446589  96.53496 


THE  YEARLY  TOTAL  ($)  =  329324.2 

THE  PUMP  ENERGY  (KWH)  =  9823.913 

THE  FAN  ENERGY  (KWH)  =  15236.21 

THE  CHILLER  ENERGY  (KWH)  =  337296.9 

THE  TOTALSYSTEM  ENERGY  (KWH)  =  362357 


58 


ATEF12.DAT  TOWER  #  12  EFFECTIVENESS  METHOD 


MONTH=  12 
TIME  START  = 
Ctrl*  7 
1=  12 

TDB(DEG  F) 

57 

GPM 

750 

CWS(DEG  F) 

80 

9.327599E-02 
PIPE  D(FT) 
.6666667 


BIN  =  57 

16:28:21  TIME  END  =  16:35:14  RUNTIME  =  6:53 
Ctr2  *  93  Ctr3  *  3 


J=  10 

TWB(DEG  F) 
52 

%GPM 

75 

CWR(DEG  F) 

84.03179 


HOURS 

24.82143 

CFM 

120000 

DELTATCOND 


TOWER 

12 

TONS 

100 

DELTATTOWER 


TOWERDUTY 


4.031794 


21.91807 


Pipevelocity (FT/SEC)  REd 
4.787051  369371.3 


pipe  frict 
1.592428E-02 


pipeL  total  PHpipe  PHcond 

PHtotal(FT  Water) 

666.6667  5.666439  5.105591 

PASSES  Condenservelocity 


PHtower 


2 

UACOND 

893100.2 

KWTONGUESS 

.9074577 

HOURS 

24.82143 

CHILLERKW 

90.74577 


2.900009 

LMTDCOND 

1.690421 

KWTONMODEL 

.9074562 

HOURSSUM 

57.5 

PUMPKW 

4.342596 


NfISENT 

.5164742 

KWHBIN 

2390.465 


FANKW 

1.218139 


11 


NFTON 

.6406521 

KWHTOT 

6136.881 


TOTALKW 

96.30651 


21.77203 

COND.TEMP 

84.4407 


THE  YEARLY  TOTAL  ($)  * 

THE  PUMP  ENERGY  (KWH)  = 

THE  FAN  ENERGY  (KWH)  = 

THE  CHILLER  ENERGY  (KWH)  = 
THE  TOTALS YSTEM  ENERGY  (KWH) 


328978.1 

9823.913 

12646.48 

337296.8 


359767.2 
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ATEF16.DAT  TOWER  #  16  EFFECTIVENESS  METHOD 


MONTH=  12  BIN  «  57 

TIME  START  =  16:36:34  TIME  END  *  16:44:47  RUNTIME  =8:13 
ctrl=  7  ctr2  =  91  ctr3  =  3 

1=  12  J=  10 

TDB(DEG  F)  TWB(DEG  F)  HOURS  TOWER 

57  52  24.82143  16 

GPM  %GPM  CFM  TONS 

750  75  120000  100 

CWS(DEG  F)  CWR(DEG  F)  DELTATCOND  DELTATTOWER  TOWERDUTY 

80  84.03179  4.031794  22.53072 

8.771837E-02 

PIPE  D(FT)  Pipevelocity (FT/SEC)  REd  pipe  frict 

.6666667  4.787051  369371.3  1.592428E-02 

pipeL  total  PHpipe  PHcond  PHtower 

PHtotal(FT  Water) 

666.6667  5.666439  5.105591  11  21.77203 

PASSES  Condenservelocity 

2  2.900009 

UACOND  LMTDCOND  NfISENT  NFTON  COND.TEMP 

893100.2  1.690421  .5164742  .6406521  84.4407 

KWTONGUESS  KWTONMODEL  KWHBIN  KWHTOT 

.9074577  .9074562  2380.363  6098.738 

HOURS  HOURS SUM 

24.82143  57.5 

CHILLERKW  PUMPKW  FANKW  TOTALKW 

90.74577  4.342596  .8111207  95.89949 


THE  YEARLY  TOTAL  ($)  =  328464.3 

THE  PUMP  ENERGY  (KWH)  =  9823.913 

THE  FAN  ENERGY  (KWH)  =  8709.041 

THE  CHILLER  ENERGY  (KWH)  =  337297.1 

THE  TOTALSYSTEM  ENERGY  (KWH)  =  355830.1 
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ATEF19.DAT  TOWER  #  19  EFFECTIVENESS  METHOD 


MONTH=  12  BIN  =  57 

TIME  START  =  16:46:15  TIME  END  »  16:55:43  RUNTIME  =9:28 
ctrl=  8  ctr2  =  104  ctr3  =  3 


1=  12 

TDB(DEG  F) 

57 

GPM 

750 

CWS(DEG  F) 

80 

8.117554E-02 
PIPE  D(FT) 
.6666667 


J=  10 

TWB(DEG  F) 
52 

%GPM 

75 

CWR(DEG  F) 

84.03179 


HOURS 

24.82143 

CFM 

120000 

DELTATCOND 


TOWER 

19 

TONS 

100 

DELTATTOWER 


TOWERDUTY 


4.031794 


23.29736 


Pipevelocity (FT/SEC)  REd 
4.787051  369371.3 


pipe  frict 
1.592428E-02 


pipeL  total  PHpipe  PHcond 

PHtotal(FT  Water) 

666.6667  5.666439  5.105591 

PASSES  Condenservelocity 


PHtower 


2 

UACOND 

893100.2 

KWTONGUESS 

.9074578 

HOURS 

24.82143 

CHILLERKW 

90.74578 


2.900009 

LMTDCOND 

1.690421 

KWTONMODEL 

.9074562 

HOURSSUM 

57.5 

PUMPKW 

4.342596 


NfISENT 

.5164742 

KWHBIN 

2379.924 


FANKW 

.7934307 


11 


NFTON 

.6406521 

KWHTOT 

6098.561 


TOTALKW 

95.88181 


21.77203 

COND.TEMP 

84.4407 


THE  YEARLY  TOTAL  ($)  = 

THE  PUMP  ENERGY  (KWH)  = 

THE  FAN  ENERGY  (KWH)  = 

THE  CHILLER  ENERGY  (KWH)  = 
THE  TOTALS Y STEM  ENERGY  (KWH) 


328486.1 

9823.913 

8849.848 

337296.8 


355970.5 
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ATFM08.DAT  CFN  80,000  JOYCE  METHOD 


MONTH=  12 

BIN  =  57 

START  TIME 

=  14:10:57 

END  TIME  =  14:51 

:45  RUNTIME  = 

=  40:48 

Ctrl=  6 

ctr2  =  80 

ctr3  *  3 

1=  12 

J=  10 

TDB(DEG  F) 

TWB(DEG 

F)  HOURS 

TOWER 

57 

52 

24.82143 

8 

GPM 

%GPM 

CFM 

TONS 

750 

75 

80000 

100 

CWS(DEG  F) 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER 

TOWERDUTY 

80 

84.03179 

4.031793 

14.52383 

.1973316 

PIPE  D(FT) 

Pipevelocity (FT/SEC)  REd 

pipe 

frict 

.6666667 

4.787051 

369371.3 

1.592428E- 

02 

pipeL  total  PHpipe 

PHcond 

PHtower 

PHtotal(FT 

Water) 

666.6667 

5.666439 

5.105591 

11 

21.77203 

PASSES 

Condenservelocity 

2 

2.900009 

UACOND 

LMTDCOND 

NfISENT 

NFTON 

COND . TEMP 

893100.2 

1.690421 

.5164742 

.6406521 

84.4407 

KWTONGUESS 

KWTONMODEL 

KWHBIN 

KWHTOT 

.9074569 

.9074562 

2386.287 

6108.143 

HOURS 

HOURSSUM 

24.82143 

57.5 

CHILLERKW 

PUMPKW 

FANKW 

TOTALKW 

90.74569 

4.342596 

1.049877 

96.13816 

THE  YEARLY  TOTAL  ($)  =  329347.2 

THE  PUMP  ENERGY  (KWH)  =  9823.913 

THE  FAN  ENERGY  (KWH)  =  7669.604 

THE  CHILLER  ENERGY  (KWH)  =  341085.2 

THE  TOTALSYSTEM  ENERGY  (KWH)  =  358578.7 


62 


ATFM10.DAT  CPM  100,000  JOYCE  METHOD 


MONTH=  12 

BIN  =  57 

START  TIME 

=  15:34:44 

END  TIME  =  16:17 

:07  RUITTIME  =  42:23 

ctrl=  6 

ctr2  *  80 

ctr3  =  3 

1=  12 

J=  10 

TDB(DEG  F) 

TWB(DEG 

F)  HOURS 

TOWER 

57 

52 

24.82143 

8 

GPM 

%GPM 

CFM 

TONS 

750 

75 

100000 

100 

CWS(DEG  F) 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER  TOWERDUTY 

80 

84.03179 

4.031794 

16.90012  .1539621 

PIPE  D(FT) 

Pipevelocity (FT/SEC)  REd 

pipe  frict 

.6666667 

4.787051 

369371.3 

1.592428E-02 

pipeL  total  PHpipe 

PHcond 

PHtower 

PHtotal(FT 

Water) 

666.6667 

5.666439 

5.105591 

11  21.77203 

PASSES 

Condenservelocity 

2 

2.900009 

UACOND 

LMTDCOND 

NfISENT 

NFTON  COND . TEMP 

893100.2 

1.690421 

.5164742 

.6406521  84.4407 

KWTONGUESS 

KWTONMODEL 

KWHBIN 

KWHTOT 

.9074571 

.9074562 

2401.928 

6167.014 

HOURS 

HOURS SUM 

• 

24.82143 

57.5 

CHILLERKW 

PUMPKW 

FANKW 

TOTALKW 

90.74571 

4.342596 

1.680003 

96.76831 

THE  YEARLY  TOTAL  ($)  =  329446 

THE  PUMP  ENERGY  (KWH)  =  9823.913 

THE  FAN  ENERGY  (KWH)  =  13997.07 

THE  CHILLER  ENERGY  (KWH)  =  338508.3 

THE  TOTALSYSTEM  ENERGY  (KWH)  =  362329.2 
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ATFM12.DAT  CFM  120,000  JOYCE  METHOD 


MONTH=  12 

BIN  =  57 

START  TIME 

=  09:38:23  END  TIME  *  10:19 

:09  RUNTIME  = 

=  40:46 

Ctrl=  7 

Ctr2  =  91 

Ctr3  *  3 

1=  12 

J=  10 

TDB(DEG  F) 

TWB(DEG  F) 

HOURS 

TOWER 

57 

52 

24.82143 

8 

GPM 

%GPM 

CFM 

TONS 

750 

75 

120000 

100 

CWS(DEG  F) 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER 

TOWERDUTY 

80 

84.03179 

4.031794 

18.92442 

.1256079 

PIPE  D(FT) 

Pipevelocity (FT/SEC)  REd 

pipe  frict 

.6666667 

4.787051 

369371.3 

1.592428E-02 

pipeL  total  PHpipe 

PHcond 

PHtower 

PHtotal(FT 

Water) 

666.6667 

5.666439 

5.105591 

11 

21.77203 

PASSES 

Condenservelocity 

2 

2.900009 

UACOND 

LMTDCOND 

NfISENT 

NFTON 

COND . TEMP 

893100.2 

1.690421 

.5164742 

.6406521 

84.4407 

KWTONGUESS 

KWTONMODEL 

KWHBIN 

KWHTOT 

.9074574 

.9074562 

2421.411 

6244.041 

HOURS 

HOURSSUM 

24.82143 

57.5 

CHILLERKW 

PUMPKW 

FANKW 

TOTALKW 

90.74574 

4.342596 

2.464895 

97.55323 

THE  YEARLY  TOTAL  ($)  =  330303.3 

THE  PUMP  ENERGY  (KWH)  =  9823.913 

THE  FAN  ENERGY  (KWH)  =  22643.94 

THE  CHILLER  ENERGY  (KWH)  =  337620.4 

THE  TOTALSYSTEM  ENERGY  (KWH)  =  370088.2 
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ATFM14.DAT 


CFM  140,000  JOYCE  METHOD 


MONTH=  12  BIN  »  57 

START  TIME  =  09:20:36  END  TIME  =  10:03:05  RUNTIME  =  42:29 
Ctrl®  7  ctr2  =  90  ctr3  =  3 


1=  12 

TDB(DEG  F) 
57 
GPM 
750 

CWS(DEG  F) 
80 

PIPE  D(FT) 
.6666667 
pipeL  total 
PHtOtal(FT 
666.6667 
PASSES 
2 

UACOND 

893100.2 

KWTONGUESS 

.9074576 

HOURS 

24.82143 

CHILLERKW 

90.74576 


J=  10 

TWB(DEG  F) 
52 

%GPM 

75 

CWR(DEG  F) 
84.03179 


2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURSSUM 
57.5 
PUMPKW 
4.342596 


HOURS 

24.82143 

CFM 

140000 

DELTATCOND 

4.031794 


NfISENT 

.5164742 

KWHBIN 

2444.862 


FANKW 

3.409683 


Pipevelocity (FT/SEC)  REd 
4.787051  369371.3 

PHpipe  PHcond 

Water) 

5.666439  5.105591 

Condenservelocity 


TOWER 

8 

TONS 

100 

DELTATTOWER  TOWERDUTY 

20.65385  .1057865 

pipe  frict 
1.592428E-02 
PHtower 


11  21.77203 


NFTON  COND.TEMP 

.6406521  84.4407 

KWHTOT 
6341.728 


TOTALKW 

98.49804 


THE  YEARLY  TOTAL  ($)  = 

THE  PUMP  ENERGY  (KWH)  = 

THE  FAN  ENERGY  (KWH)  = 

THE  CHILLER  ENERGY  (KWH)  = 
THE  TOTALS YSTEM  ENERGY  (KWH) 


331787.7 

9823.913 

34027.18 

337324.1 


381175.2 
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ATFM16.DAT  CFM  160,000  JOYCE  METHOD 


MONTH=  12 

BIN  =  57 

START  TIME  = 

10:04:59  END 

TIME  *  10:48 

:51  RUNTIME  =  43:52 

ctrl=  8 

ctr2  =  100 

ctr3  *  3 

1=  12 

J=  10 

TDB(DEG  F) 

TWB(DEG  F) 

HOURS 

TOWER 

57 

52 

24.82143 

8 

GPM 

%GPM 

CFM 

TONS 

750 

75 

160000 

100 

CWS(DEG  F) 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER  TOWERDUTY 

80 

84.03179 

4.031794  22.13515 

9.127153E-02 

PIPE  D(FT) 

Pipevelocity (FT/SEC)  REd 

pipe  frict 

.6666667 

4.787051 

369371.3 

1.592428E-02 

pipeL  total 

PHpipe 

PHcond 

PHtower 

PHtotal(FT  Water) 

666.6667 

5.666439 

5.105591 

11  21.77203 

PASSES 

Condenservelocity 

2 

2.900009 

UACOND 

LMTDCOND 

NfISENT 

NFTON  COND.TEMP 

893100.2 

1.690421 

.5164742 

.6406521  84.4407 

KWTONGUESS 

KWTONMODEL 

KWHBIN 

KWHTOT 

.9074579 

.9074562 

2472.463 

6462.413 

HOURS 

HOURSSUM 

24.82143 

57.5 

CHILLERKW 

PUMPKW 

FANKW 

TOTALKW 

90.74579 

4.342596 

4.521638 

99.61002 

THE 

YEARLY  TOTAL  ($)  = 

333746.1 

THE 

PUMP  ENERGY  (KWH)  = 

9823.913 

THE 

FAN  ENERGY  (KWH)  = 

48544.85 

THE 

CHILLER  ENERGY  (KWH)  = 

337296.6 

THE 

TOTALS Y STEM  ENERGY  (KWH) 

= 

395665.4 
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ATFM18.DAT  CFM  180,000  JOYCE  METHOD 


MONTH=  12 

BIN  =  57 

START  TIME 

=  13:28:34  END 

TIME  =  14:09 

:48  RUNTIME  = 

41:14 

Ctrl=  8 

ctr2  =  94 

ctr3  =  3 

1=  12 

J=  10 

TDB(DEG  F) 

TWB(DEG  F) 

HOURS 

TOWER 

57 

52 

24.82143 

8 

GPM 

%GPM 

CFM 

TONS 

750 

75 

180000 

100 

CWS(DEG  F) 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER 

towerditty 

80 

84.03179 

4.031794 

23.40676 

.08P2768 

PIPE  D(FT) 

Pipevelocity (FT/SEC)  REd 

pipe  frict 

.6666667 

4.787051 

369371.3 

1.592428E-02 

pipeL  total  PHpipe 

PHcond 

PHtower 

PHtotal(FT 

Water) 

666.6667 

5.666439 

5.105591 

11 

21.77203 

PASSES 

Condenserve loc i ty 

2 

2.900009 

UACOND 

LMTDCOND 

NfISENT 

NFTON 

COND . TEMP 

893100.2 

1.690421 

.5164742 

.6406521 

84.4407 

KWTONGUESS 

KWTONMODEL 

KWHBIN 

KWHTOT 

.9074581 

.9074562 

2504.454 

6608.818 

HOURS 

HOURSSUM 

24.82143 

57.5 

CHILLERKW 

PUMPKW 

FANKW 

TOTALKW 

90.7458 

4.342596 

5.810467 

100.8989 

THE  YEARLY  TOTAL  ($)  =  336225.5 

THE  PUMP  ENERGY  (KWH)  =  9823.913 

THE  FAN  ENERGY  (KWH)  =  66800.08 

THE  CHILLER  ENERGY  (KWH)  =  337297.3 

THE  TOTALSYSTEM  ENERGY  (KWH)  =  413921.3 
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AFMEF08.DAT  CFM  80,000  EFFECTIVENESS  METHOD 


MONTH=  12 

BIN  =  57 

TIME  START  = 

14:16:25  TIME  END  *  14:24 

:02  RUNTIME  =  7:37 

ctrl=  6 

ctr2  s=  85 

ctr3  *  3 

1=  12 

J=  10 

TDB(DEG  F) 

TWB(DEG  F) 

HOURS 

TOWER 

57 

52 

24.82143 

8 

GPM 

%GPM 

CFM 

TONS 

750 

75 

80000 

100 

CWS(DEG  F) 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER  TOWERDUTY 

80 

84.03179 

4.031794 

15.58475  .1763346 

PIPE  D(FT) 

Pipevelocity (FT/SEC)  REd 

pipe  frict 

.6666667 

4.787051 

369371.3 

1.592428E-02 

pipeL  total 

PHplpe 

PHcond 

PHtower 

PHtotal(FT  Water) 

666.6667 

5.666439 

5.105591 

11  21.77203 

PASSES 

Condenservelocity 

2 

2.900009 

UACOND 

LMTDCOND 

NfISENT 

NFTON  COND.TEMP 

893100.2 

1.690421 

.5164742 

.6406521  84.4407 

KWTONGUESS 

KWTONMODEL 

KWHBIN 

KWHTOT 

.907457 

.9074562 

2383.514 

6099.48 

HOURS 

HOURS SUM 

• 

24.82143 

57.5 

CHILLERKW 

PUMPKW 

FANKW 

TOTALKW 

90.7457 

4.342596 

.9381654 

96.02646 

THE  YEARLY  TOTAL  ($)  =  328801.4 

THE  PUMP  ENERGY  (KWH)  =  9823.913 

THE  FAN  ENERGY  (KWH)  *  7206.588 

THE  CHILLER  ENERGY  (KWH)  =  339184.5 

THE  TOTALS YSTEM  ENERGY  (KWH)  =  356215 
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AFMEF10.DAT  CFM  100,000  EFFECTIVENESS  METHOD 


MONTH=  12  BIN  =57 

TIME  START  =  14:28:05  TIME  END  =  14:36:57  RUNTIME  =  8:52 


Ctrl=  6 
1=  12 

TDB(DEG  F) 
57 
GPM 
750 

CWS(DEG  F) 
80 

PIPE  D(FT) 
.6666667 
pipeL  total 
PHtotal(FT 
666.6667 
PASSES 
2 

UACOND 

893100.2 

KWTONGUrSS 

.9074572 

HOURS 

24.82143 

CHILLERKW 

90.74572 


Ctr2  =  81 

J=  10 

TWB(DEG  F) 
52 

%GPM 

75 

CWR(DEG  F) 
84.03179 


2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURSSUM 
57.5 
PUMPKW 
4.342596 


Ctr3  =  3 

HOURS 

24.82143 

CFM 

100000 

DELTATCOND 

4.031794 


NfISENT 

.5164742 

KWHBIN 

2396.03 


FANKW 

1.442381 


Pipevelocity (FT/SEC)  REd 
4.787051  369371.3 

PHpipe  PHcond 

Water) 

5.666439  5.105591 

Condenservelocity 


TOWER 

8 

TONS 

100 

DELTATTOWER  TOWERDUTY 

18.4128  .1321855 

pipe  frict 
1.592428E-02 
PHtower 


11 


NFTON 

.6406521 

KWHTOT 

6148.765 


TOTALKW 

96.53069 


21.77203 

COND.TEMP 

84.4407 


THE  YEARLY  TOTAL  ($)  = 

THE  PUMT"  ENERGY  (KWH)  = 

THE  FAN  ENERGY  (KWH)  = 

THE  CHILLER  ENERGY  (KWH)  = 
THE  TOTALSYSTEM  ENERGY  (KWH) 


329002.8 

9823.913 

12707.91 

337642.3 


360174.1 
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AFMEF12.DAT  CFM  120/000  EFFECTIVENESS  METHOD 


MONTH=  12 

BIN  =  57 

TIME  START  = 

12:46:37  TIME  END  =  12:53: 

25  RUNTIME  *  6:48 

ctrl=  7 

Ctr2  =  92 

Ctr3  *  3 

1=  12 

J=  10 

TDB(DEG  F) 

TWB(DEG  F) 

HOURS 

TOWER 

57 

52 

24.82143 

8 

GPM 

%GPM 

CFM 

TONS 

750 

75 

120000 

100 

CWS(DEG  F) 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER  TOWERDUTY 

80 

84.03179 

4.031794 

20.85859  .1036575 

PIPE  D(FT) 

Pipevelocity (FT/SEC)  REd 

pipe  frict 

.6666667 

4.787051 

369371.3 

1.592428E-02 

pipeL  total 

PHpipe 

PHcond 

PHtower 

PHtotal(FT  Water) 

666.6667 

5.666439 

5.105591 

11  21.77203 

PASSES 

Condenservelocity 

2 

2.900009 

UACOND 

LMTDCOND 

NfISENT 

NFTON  COND . TEMP 

893100.2 

1.690421 

.5164742 

.6406521  84.4407 

KWTONGUESS 

KWTONMODEL 

KWHBIN 

KWHTOT 

.9074575 

.9074562 

2410.719 

6211.164 

HOURS 

HOURSSUM 

24.82143 

57.5 

CHILLERKW 

PUMPKW 

FANKW 

TOTALKW 

90.74575 

4.342596 

2.034147 

97.12249 

THE  YEARLY  TOTAL  ($)  =  329947.8 

THE  PUMP  ENERGY  (KWH)  =  9823.913 

THE  FAN  ENERGY  (KWH)  =  20089.43 

THE  CHILLER  ENERGY  (KWH)  =  337312.3 

THE  TOTALSYSTEM  ENERGY  (KWH)  =  367225.7 
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AFMEF14.DAT  CFM  140,000  EFFECTIVENESS  METHOD 


MONTH=  12  BIN  »  57 

TIME  START  =  12:54:56  TIME  END  =  13:03:20  RUNTIME  =  8:24 
ctrl=  8  ctr2  =  104  ctr3  *  3 


1=  12 

TDB(DEG  F) 

57 

GPM 

750 

CWS(DEG  F) 

80 

8.405238E-02 
PIPE  D(FT) 
.6666667 


J=  10 

TWB(DEG  F) 
52 

%GPM 

75 

CWR(DEG  F) 

84.03179 


HOURS 

24.82143 

CFM 

140000 

DELTATCOND 


TOWER 

8 

TONS 

100 

DELTATTOWER 


TOWERDUTY 


4.031794 


22.95394 


Pipevelocity (FT/SEC)  REd 
4.787051  369371.3 


pipe  frict 
1.592428E-02 


pipeL  total  PHpipe 
PHtotal(FT  Water) 


PKcond 


666.6667 

PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074579 

HOURS 

24.82143 

CHILLERKW 

90.74579 


5.666439  5.105591 

Condenservelocity 
2.900009 


LMTDCOND 

1.690421 

KWTONMODEL 

.9074562 

HOURSSUM 

57.5 

PUMPKW 

4.342596 


NfISENT 

.5164742 

KWHBIN 

2427.475 


FANKW 

2.709155 


PHtower 


11  21.77203 


NFTON  COND . TEMP 

.6406521  84.4407 

KWHTOT 
6288.67 


TOTALKW 

97.79754 


THE  YEARLY  TOTAL  ($)  = 

THE  PUMP  ENERGY  (KWH)  = 

THE  FAN  ENERGY  (KWH)  = 

THE  CHILLER  ENERGY  (KWH)  = 
THE  TOTALSYSTEM  ENERGY  (KWH) 


331261.7 

9823.913 

29775.94 

337297 


376896.8 
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AFMEF16.0AT  CFM  160,000  EFFECTIVENESS  METHOD 


MONTH=  12 

BIN  =  57 

TIME  START 

=  11:06:16  TIME  END  =  11:13 

:18  RUNTIME  =7:02 

ctrl=  8 

Ctr2  =  105 

Ctr3  **  3 

1=  12 

J=  10 

TDB(DEG  F) 

TWB(DEG  F) 

HOURS 

TOWER 

57 

52 

24.82143 

8 

GPM 

%GPM 

CFM 

TONS 

750 

75 

160000 

100 

CWS(DEG  F) 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER  TOWERDUTY 

80 

84.03179 

4.031794 

24.73064  .0700314 

PIPE  D(FT) 

Pipevelocity (FT/SEC)  REd 

pipe  frict 

.6666667 

4.787051 

369371.3 

1.592428E-02 

plpeL  total  PHpipe 

PHcond 

PHtower 

PHtotal(FT 

Water) 

666.6667 

5.666439 

5.105591 

11  21.77203 

PASSES 

Condenservelocity 

2 

2.900009 

UACOND 

LMTDCOND 

NfISENT 

NFTON  COND.TEMP 

893100.2 

1.690421 

.5164742 

.6406521  84.4407 

KWTONGUESS 

KWTONMODEL 

KWHBIN 

KWHTOT 

.9074581 

.9074562 

2446.345 

6383.272 

HOURS 

HOURSSUM 

24.82143 

57.5 

CHILLERKW 

PUMPKW 

FANKW 

TOTALKW 

90.7458 

4.342596 

3.469392 

98.55779 

THE  YEARLY  TOTAL  ($)  =  332982.5 

THE  PUMP  ENERGY  (KWH)  =  9823.913 

THE  FAN  ENERGY  (KWH)  =  42322.18 

THE  CHILLER  ENERGY  (KWH)  =  337297.3 

THE  TOTALSYSTEM  ENERGY  (KWH)  =  389443.3 
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AFMEF18.DAT  CFM  180«000  EFFECTIVENESS  METHOD 

MONTH=  12 

BIN  =  57 

TIME  START  = 

12:18:54  TIME  END  =  12:27: 

39  RUNTIME  =8:45 

ctrl=  9 

ctr2  =  116 

Ctr3  *  3 

1=  12 

J=  10 

TDB(DEG  F) 

TWB(DEG  F) 

HOURS 

TOWER 

57 

52 

24.82143 

8 

GPM 

%GPM 

CFM 

TONS 

750 

75 

180000 

100 

CWS(DEG  F) 

CWR(DEG  F) 

DELTATCOND 

DELTATTOWER  TOWERDUTY 

80 

84.03179 

4.031795  26.22048 

5.973898E-02 

PIPE  D(FT) 

Pipevelocity (FT/SEC)  REd 

pipe  frict 

.6666667 

4.787051 

369371.3 

1.592428E-02 

pipeL  total 

PHpipe 

PHcond 

PHtower 

PHtotal(FT  Water) 

666.6667 

5.666439 

5.105591 

11  21.77203 

PASSES 

Condenserve loc i ty 

2 

2.900009 

UACOND 

LMTDCOND 

NfISENT 

NFTON  COND.TEMP 

893100.2 

1.690421 

.5164742 

.6406521  84.4407 

KWTONGUESS 

KWTONMODEL 

KWHBIN 

KWHTOT 

.9074584 

.9074562 

2467.557 

6497.896 

HOURS 

HOURSSUM 

24.82143 

57.5 

CHILLERKW 

PUMPKW 

FANKW 

TOTALKW 

90.74583 

4.342596 

4.323931 

99.41236 

THE  YEARLY  TOTAL  ($)  =  335188.1 

THE  PUMP  ENERGY  (KWH)  *  9823.913 

THE  FAN  ENERGY  (KWH)  =  58278.78 

THE  CHILLER  ENERGY  (KWH)  =  337297.4 

THE  TOTALSYSTEM  ENERGY  (KWH)  =  405400.2 
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